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SUMMARY 


General analytical methods for the prediction of the attitude motion of a gravity 
gradient stabilized spacecraft are developed for an elliptic orbit. Both planar and 
three-dimensional motion are treated. 

The linearized small eccentricity equation describing the single degree of 
freedom in-plane motion (pitch motion) is a Mathieu equation. Combinations of 
spacecraft inertia parameters and orbit eccentricities for which the spacecraft 
pitch motion is unbounded are plotted from tables of the Mathieu function. Three 
regions of unstable motion are found. 

The linearized planar equation for larger eccentricities is treated by means 
of the asymptotic expansion methods (in powers of the eccentricity) of Krylov, 
Bogoliubov, andMitropolsky. A more accurate determination of the three regions 
of unstable motion found from Mathieu’s equation is obtained. 

Finally an analysis based on the Hamiltonian formulation of the equations of 
motion demonstrates the effect of both nonlinear and parametric resonances on the 
three degree of freedom motion of the spacecraft. Canonical transformations and 
the method of averaging are used to determine combinations of eigenfrequencies of 
the normal modes which lead to various types of resonances. This is an extension 
of an earlier work by Breakwell and Pringle in which six resonant frequency com- 
binations were found. Twenty- six additional resonance combinations are found of 
which sixteen lead to unbounded motion in one or more of the modes and six lead to 
significant interchanges of energy between various modes. 

It is concluded first that the attitude motion of a spacecraft in an eccentric 
orbit has characteristics significantly different from those of the corresponding 
motion in a circular orbit and second that theoretical methods are adequate to pre- 
dict, and hence avoid, the occurrence of unstable modes during the design of gravity 
gradient spacecraft. 
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ON THE LIBRATION OF A GRAVITY GRADIENT 
STABILIZED SPACECRAFT IN AN ECCENTRIC ORBIT 

by 

Albert J. Fleig, Jr. 


CHAPTER 1 
INTRODUCTION 


Historical Background 

Analytic investigations of the effect of an inverse-square gravitational field on a triaxial sat- 
ellite have been published over a span of nearly 300 years. In 1686, Newton (Reference 1) wrote: 

", . . the figure of the moon would be a spheroid, whose greatest diameter produced would pass 
through the center of the earth. . . . Hence it is that the same face of the moon always is turned 
toward the earth; nor can the body of the moon possibly rest in any other position, but would re- 
turn always by a libratory motion to this situation. . . According to Routh (Reference 2), during 
the next 200 years the attention of D'Alembert, Lagrange, Laplace, and Poisson was directed to- 
ward the gravity- gradient effect on the motion of the moon. By the beginning of the 20th century, 
correct equations for this torque were routinely published as in Plummer (Reference 3). One of 
the goals of these early investigators was "selenodesy," i.e., making inferences about the struc- 
ture of the moon from observation of its motion. Present emphasis, apparently starting with a 
proposal for a hinged satellite by Breakwell and Roberson in 1954, centers on the other half of this 
problem, predicting spacecraft motion from a knowledge of its structure. 

The particular problem that forms the subject of this investigation— the motion of a gravity- 
gradient- stabilized spacecraft in an eccentric orbit— was first discussed by Baker (Reference 4). 

He commented that motion about the orbit normal (i.e. pitch) could be approximately represented by 
a linearized equation that was similar to a standard Mathieu equation with a forcing term added. 
Schrello (Reference 5) in his investigations of the gravity-gradient disturbance of aerodynamically 
stable spacecraft formulated the pitch motion problem in terms of Hill's equation. DeBra (Refer- 
ence 6) published a numerical investigation of the effect of orbit eccentricity on attitude stability. 

A number of Russian authors, particularly Beletskii (References 6, 8, and 9) and Chernous'ko (Ref- 
erence 10), made an analysis in which they applied asymptotic methods to the problem. Kane (Ref- 
erence 11) applied Floquet theory to a somewhat restricted three-dimensional analysis and discovered 


*In its original form, this paper was submitted to the Catholic University of America in partial fulfillment of the requirements for the 
degree of Doctor of Philosophy. 
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apparent instabilities in the nominally stable range of spacecraft parameters. Breakwell and 
Pringle (Reference 12) explained this phenomenon in terms of a nonlinear resonance between the 
pitch and yaw frequencies. 

This dissertation contains several analytical approaches to the problem of determining the 
librational motion of rigid spacecraft in an eccentric orbit about an ideal spherical earth. The 
analysis is primarily concerned with the effect of orbit eccentricity on spacecraft motion; no other 
perturbations are considered. The orbital motion (motion of the spacecraft’s center of mass) is 
considered known a priori. 


Objectives 

The exact equations of motion for a gravity- gradient- stabilized spacecraft in an eccentric orbit 
are coupled* time-varying* nonlinear, and insoluble in closed form. Thus analytic efforts have been 
concentrated on a number of special cases such as the linearized equations for three-dimensional 
motion in a circular orbit, the nonlinear equation for pitch motion in a circular orbit, the linear 
equation for pitch motion in an eccentric orbit, etc. Since a complete analytical solution is unavail- 
able, parameter selection for any particular spacecraft is based on a combination of linear analysis 
and numerical simulation. Yet linear analysis is unable to reveal many significant features of the 
equations of motion; furthermore, the basic phenomena underlying the results obtained from a nu- 
merical simulation are seldom explained. 

The major goal of this dissertation is to provide an explicit analytical basis for the effect of 
orbit eccentricity on spacecraft attitude. Two separate approaches to the problem are made. First, 
the asymptotic methods of Krylov, Bogoliubov, and Mitropolsky are applied to the planar-pitch 
problem to obtain stability regions and, in the stable regions, to determine the forced responses. 
Then the nonlinear, three-dimensional problem is discussed, using canonical transformation and 
the method of averaging to obtain a locus of spacecraft inertia parameters that result in resonant 
motion. 

Other objectives are to provide increased understanding of one of the mechanisms by which 
nonlinear terms can produce unexpected results in problems with many degrees of freedom, and a 
basis for application of the two approaches to particular mechanizations of the basic gravity- 
gradient- stabilization concept (e.g., to consider specific damping concepts). 

Outline 

Chapter 2 discusses the linearized equation of spacecraft motion about the pitch axis (i.e. orbit 
normal), with the expansions of orbital parameters restricted to first order in eccentricity, e. This 
is the simplest form of the problem. There are several possible variations involving choices for 
reference axis and independent variable. Each of the choices leads to a Mathieu-type equation with 
a forcing term, but the final equations are not identical. Regions of stable motion are plotted as 
functions of orbit eccentricity and spacecraft inertia ratio for each of the formulations. The stability 
plots are similar if e is very small but diverge rapidly as e increases. 
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Chapter 3 develops an asymptotic solution for the in-plane motion with a more accurate ex- 
pansion for the orbital elements. When the frequency of oscillation about the pitch axis is 1/2, 1, 
or 3/2 cycles per orbit, there is a parametric excitation of the pitch motion, as in Chapter 2. The 
basic asymptotic solution for this problem determines the magnitude of the induced oscillation 
about the nominal equilibrium axis (i.e., local vertical). However, this form of asymptotic solution 
fails in the vicinity of each parametric resonance. A second type of solution is developed for each 
of the resonant frequencies. These solutions are used to obtain boundaries for stable motion sim- 
ilar to those developed in Chapter 2. This completes the analytical treatment of the pitch motion 
case with one degree of freedom. The next two chapters are concerned with the nonlinear equa- 
tions for the full three degrees of freedom. 

A different approach, first applied to this problem by Breakwell and Pringle (Reference 12) is 
used in the discussion of spacecraft motion with three degrees of freedom. A linearized description 
of the motion in terms of normal modes leads to three uncoupled simple -harmonic oscillators in 
normal coordinates. Nonlinear terms and time-varying terms due to orbital eccentricity tend to 
cause changes in the fundamental uncoupled oscillations. These changes occur slowly in compari- 
son with the time for a single cycle of any of the normal modes, and only become significant when 
there is a resonant relationship between the frequency of the forcing terms and the appropriate 
normal mode(s). 

Chapter 4 develops the Hamiltonian for the system to the third order in generalized coordin- 
ates and momenta and with orbital motion represented in expansions to fourth order in e; also, a 
series of cannonical transformations that first decouple the roll-yaw equations and then introduce 
cyclic coordinates. The resulting equations of motion represent a set of three harmonic oscillators 
perturbed by second-order terms and forcing functions. 

When the method of averaging is applied to the Hamiltonian in Chapter 5, the effect of the non- 
linear coupling terms and the forcing terms becomes apparent. The amplitude and phase of the 
fundamental oscillations vary only slightly from their initial values except when there are reson- 
ances between the perturbing terms and the fundamental modes. Approximately 30 resonant rela- 
tions exist (when the Hamiltonian expansion is limited as noted above); the locus of each of the 
resonances is plotted in this chapter. It is possible to define regions about each of the resonances 
in which the motion would change appreciably. The conditions which define these regions are given 
for each of the possible resonances. 

The appendices contain a discussion of the coupling between orbital motion and attitude motion 
and the specific extensions of the asymptotic formulas used in Chapter 3 for the third-order 
solutions. 

Comments on Linearization and Stability 

Linear differential equations are so convenient, and dynamics problems starting with the pen- 
dulum are so often linearized, that linearization is frequently resorted to without adequate consid- 
eration of its validity. Speaking somewhat broadly, it is valid to linearize equations if the resulting 
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solution does not differ appreciably from the solution to the exact equations. The validity is fre- 
quently determined by intuition, or comparison with other similar equations; however, a more cor- 
rect justification can be obtained by comparing the analytic results for the linear equation with 
those obtained from integration of the exact equations performed, for instance, on an analog or dig- 
ital computer. There are rigorous mathematical conditions under which linearization is valid and 
these are discussed by Struble (Reference 13) and Cesari (Reference 14) among others. Unfortu- 
nately, these are only sufficient conditions for linearization to succeed or fail, not necessary ones. 

In general, if the equilibrium solution to the linear portion of the equation is asymptotically stable 
and the nonlinear portion meets certain conditions, the equilibrium solution of the complete equa- 
tion is also asymptotically stable. When the solution to the linear equation is stable but not asymp- 
totically so (as in this problem), no information concerning the adequacy of the linear solution is 
obtained. In fact, in the last chapter of this dissertation a number of unstable regions are found 
for the nonlinear equations in regions that satisfy all the stability criteria for the linearized equations. 

According to Stoker (Reference 15), stability is a concept that means different things to dif- 
ferent people. In this dissertation a linear system will be defined as asymptotically stable if all 
the characteristic numbers of the Floquet form of solution (Struble, Reference 13) have negative 
real parts, stable if no root has a positive real part, and unstable if any root has a positive real 
part. This is contrary to the standard servomechanism practice, in that an undamped bounded 
oscillator would be stable under the above definition. In discussing stability for nonlinear systems 
the concept of Liapunov stability (Reference 14) is usually satisfactory. Note that a conclusion of 
instability based on approximate equations guarantees only that the motion will increase until higher- 
order terms not included in the analysis become important. 

In formulating approximate equations there is no clear-cut rule as to the number or order of 
terms that should be retained. Similarly with asymptotic solutions there is a question as to the 
order of approximation to which the solution should be continued. The inclusion of additional terms 
in the equations or solutions will usually have one of two effects. Either the extension will intro- 
duce new phenomena into the solution, (e.g., including the nonlinear term in van-der Pol T s equation), 
or it will increase the precision of the solution of a particular problem (e.g., continuing the ex- 
pansion of sin x for a simple pendulum). In this dissertation the intent is to include any reasonable 
extension that effects a basic change in the solution, but not to add terms solely to increase the 
precision of a result. 

Notation and Definitions 

Four right-handed orthogonal coordinate frames shown in Figure 2-2, are used in this disserta- 
tion as follows: 

[E] = Ej, E 2 , E 3 An inertial system with origin at the earth’s center with: 

E x , a unit vector towards perigee of the spacecraft orbit, 

E 3 , a unit vector parallel to the spacecraft- orbit angular-momentum 
vector, and 

^2 = ^3 X ^ 1 * 
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[£] = 5i.Ca.Cs 


[S'] 




(e) - e r e 2 , e 3 


A local vertical system with origin at the spacecraft center of mass, 
with: 

X l9 a unit vector parallel to a vector from the earth center of mass 
to the spacecraft center of mass, 

X 39 a unit vector parallel to e 3 , and 

C 2 = ^3 >< V 

A mean vertical system with origin at the spacecraft center of mass, 
with: 

Xi 9 a unit vector parallel to a vector from the earth center of mass to 
that of an imaginary spacecraft in a circular orbit with the same 
period and the same phase as the subject spacecraft, 

X 3 ' , a unit vector parallel to E 3 , and 

A principal body-axis system, with: 
e x , a unit vector along the axis of least inertia, 
e 2 , a unit vector along the intermediate axis, and 
e 3 , a unit vector along the axis of maximum inertia. 


Notation relating to coordinate transformations: 

[A] = direction cosine matrix 

A 1 (77) = direction cosine matrix for a rotation through an angle 77 about the i th axis. 


The coordinate frames are located by the following rotation matrices: 


Local vertical— inertial 


[Cl 3 A 3 (v) [E] 


where v is the true anomaly and 


A 3 (v) 


cos v sin v 0 

~ sin v cos v 0 

0 0 1 


Mean vertical— inertial 


[£'] = A 3 (M)E 


where M is the mean anomaly. 
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Local vertical— mean vertical 


[£] = a 3 (v-m)U'] 


Body— local vertical 


[?] - A 1 ( 0 )A 2 ( 0 )A 3 U 3 = [A] [£] 


where <t>, t>, and \jj are Euler angles designated as pitch, roll, and yaw respectively and 


[A] 


cos <t cos 6 

cos (56 sin 6 sin 0 - sin 0 cos 0 
cos 0 sin # cos 0 + sin 0 sin 0 


sin 0 cos 6 

s in 0 s in 6 sin 0 + cos 0 cos 0 
sin 0sin 6 cos 0- cos 0 simp 


- sin 0 
cos 0 sin 0 
cos 6 COS 0 


Body— mean vertical 


[c] " A 1 ( 0 )A 2 (6?) A 3 ( 0 ') U] 


where <j> ~ d>' + v -M. 

Conventions and Units 

The unit of time is chosen so that the mean angular motion, n, is one, thus M = t. 

Units of mass and length are chosen such that i 3 = 1; and dimensionless parameters r : and 
r 2 are defined as 

i x i 2 

r , s r 3 ’ r 2 s r 3 ■ 

A superscript dot will indicate a derivative taken with respect to the independent variable 
which may be either t or v. 

A subscript zero indicates an initial value. 
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CHAPTER 2 


PLANAR PITCH LIBRATION AS A MATHIEU FUNCTION 


Introduction 

A complete formulation of the equations of motion for a rigid body in an inverse- square force 
field yields six coupled, nonlinear, second-order differential equations. The motion of a gravity- 
gradient- stabilized spacecraft cannot be determined analytically when presented in this exact form. 
However, there are possible approximations that lead to significant results. This chapter intro- 
duces a series of assumptions that reduce the problem to consideration of the small-angle pitch 
librations of a spacecraft moving in a known orbit of small eccentricity, and from there to a form 
of Mathieu T s equation. Even this highly restrictive formulation yields interesting results not ob- 
tained from circular-orbit analysis. 

Beletskii (Reference 8) and Baker (Reference 4) noted the connection between an approximate 
form of the in-plane pitch- libration problem defined below and a Mathieu equation. Beletskii dis- 
cussed an equation formulated with true anomaly as the independent variable, and Baker expressed 
the problem with time as the independent variable. This chapter demonstrates the differences that 
result from these two choices. A set of graphs are developed herein which present the stability 
boundaries associated with each of the approaches as a function of the inertia ratios of the space- 
craft and the orbital eccentricity. Finally the effect of velocity-dependent damping on spacecraft 
stability is shown. 

The first assumption for this and all the following chapters is that the motion of the space- 
craft’s center of mass (orbital motion) is not affected by motion of the spacecraft about its center 
of mass (attitude motion). With this widely accepted simplification, the equations representing 
orbital motion are identical with those for a point mass in an inverse- square force field. 

The converse of this first assumption is not true, and it is necessary to include the orbital 
parameters as known functions of time when solving for the spacecraft attitude motion. In essence, 
this is equivalent to assuming that energy can be transferred from orbital motion to attitude motion 
without affecting the orbital motion. Actually, as shown by Beletskii (Reference 9), the relative 
magnitudes of the appropriate terms make this a reasonable approximation. 

The next assumption is that the equations can be linearized when applied to small-angle motion. 
The linearized equations for roll and yaw do not involve terms in the pitch motion, and vice versa; 
pitch motion may be considered independent of roll-yaw motion. 

The coupled roll-yaw equations are homogeneous for both circular and eccentric orbits, and 
there is a solution with 0 for either case. This chapter and Chapter 3 are concerned 
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with the equation describing pitch motion under the assumption that the roll and yaw variables are 
identically zero. Thus the exact problem is reduced to consideration of the single linear second- 
order differential equation for pitch motion, with the orbital parameters known functions of time. 

In this chapter the equation for pitch libration will be formulated in several different ways. 

The variations originate from the different choices of reference frame and independent variable, 
even when the problem is restricted, as discussed above. Each of the formulations is limited to 
exclude terms higher than first order in eccentricity, and the resulting equations are transformed 
into the canonical form of Mathieu's equation. The nature of solutions to Mathieu f s equation is well 
known and regions of bounded and unbounded motion for the pitch libration are determined from 
these known solutions. 

Chapter 3 presents a somewhat more general approach to the same pitch libration problem. 

The equation is limited to exclude terms higher than third order in eccentricity, and the asymptotic 
methods of Krylov, Bogoliubov, and Mitropolsky are used to obtain regions of bounded and un- 
bounded motion. 


Equations of Motion 

The linearized equations of motion for a rigid, gravity -stabilized spacecraft moving in an ec- 
centric orbit are well known. In the notation of this dissertation they are: 


and 


r 1 0+(l-r 2 )v 2 0+(l-r 2 -r 1 )v6 J - 0, 

•V 7 + ( 1_ r i)( c ' 2 + 7t) e ‘ ( 1 " r 2~ r 1 )^ " 0 


When a circular orbit is assumed, these equations become: 


r i^ + ( 1 _ r 2 ) ^ + ( 1 " r 2“ r i) f " = 0 ' 

r 2 6 J + 4(l-r 1 )(9-(l-r 2 -r 1 )i/' ' 0 , 


and 


i + 3 ( r 2 - r 1 ) ^ - 0. 

One solution of this set of equations is 0 = r 0 = 0, and this solution can be shown to be a stable 
one for any values of r and r 2 in the "Lagrange" or TT Delp tT regions as shown in Figure 2-1 (after 
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Debra). Pringle (Reference 16) has demonstrated, by application of Liapunov’s second method, 
that spacecraft in the Lagrange region are stable only if there is any form of pervasive energy dis- 
sipation in the system. Therefore there are no practical spacecraft whose parameters fall in the 
Delp region, and this region will not be considered in the remainder of the dissertation. 

In much the same manner, 0 = \p = 0 is a solution of the two linearized roll-yaw equations for 
motion in an elliptic orbit. It is true that this solution cannot be proven stable throughout the entire 
Lagrange region. Nevertheless, a number of authors (e.g., Baker (Reference 4), Schrello (Ref- 
erence 5), Beletskii (Reference 8), Moran (Reference 17), Schecter (Reference 18), etc.) have con- 
sidered variations of the problem in which roll and yaw are assumed zero and only pitch motion is 
considered. This class of motion is commonly referred to as the T ’in-plane (or planar)" pitch- 
motion problem. 

Chapter 5 shows that for some portions of the Lagrange regions it is incorrect to assume: 

(1) that the equations can be linearized for small-angle motion, or (2) that 6 = 0 = 0 is a stable 
solution to the roll-yaw equations. However, a number of interesting points can be observed by 
investigation of pitch subject to these restrictions; both this and the next chapter are devoted to 
such an analysis. 

The previous equations are obtained in Appendix A by starting with an exact formulation of the 
problem and introducing simplifications into the equations. A different approach to the same result 
would be to write exact equations for a simplified problem; this would help in visualizing the effect 
of a change in reference systems from the local-vertical to the mean-vertical system. 

Consider a spacecraft moving in an eccentric orbit (Figure 2-2) and subject to a torque r given 

by 


T 





or, since 


therefore 


<£' - + v - M , 


T 



r 2 )(V- (V-M)] . 


The angular momentum of the spacecraft about its center of mass is either I 3 (<£+ v) or l 3 (^' + 1WI) 
depending on the choice of reference frame. Thus when the rate of change of angular momentum is 
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equated to the applied torque for motion referenced to the local-vertical frame, the result is 


3/^6 

* + 7? ( r 2 ' r l) <t> = ~ V • 


(2-1) 


Similarly, when the motion is referred to the mean-vertical frame the result is 


+ ( r 2 - r i) {*' ~ ( v -M)] ° 0 ; (2-2) 

this is the equation studied by Baker (Reference 4). 

Both Equations 2-1 and 2-2 are written with time as the independent variable. An alternative 
formulation in which the independent variable is changed to the true anomaly is of interest. The 
change of independent variables is accomplished by noting that 


/dv\* 

d 4. 

d 

/dv ■ 

\ dv 

\dt / 

+ dv 

dv 

Vdt t 

/dt 


and from celestial mechanics (see Appendix A) 


/dv\ 2 _ 

\dt / ~ — 3 (1 +e cos V) 


and 


d / dv \ dv _ 
dv \dt / dt " ~ “T 


Thus Equation 2-1 becomes, after some simplification, 


d 2 d 2e s in v dd ^( r 2 r i) _ 2e sin v 

c j v 2 1 4 e cos v dv + 1 + e cos v ~ 1 + e cos v 


(2-3) 


There is a standard transformation that removes the first derivative term in an equation of this 
form. Let 


<J>(v) - d(v) ( 1 + e cos v) . 
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Then 


V -*T' 1 _ n 

dv* " dv ( 1 + e cos v )~ + Oe sin v ( 1 + e cos v) 1 , 

d 2 <j> _ d 2 0 , dO 

^ 2 - * ( 1 + e cos v) 1 + 2 ^ esin v(l + e cos v) £ 

+ <D e cos v(l + e cosv)' 2 + <I>2e 2 sin 2 v(l+e cos v)~ 3 

and Equation 2-3 becomes 


d 2 <t> 3 (r 2 -rj +e cosv _ 

^2 + 1+ecosT ® - 2e sin v . (2-4) 

This is the form of the equation that has been studied by Beletskii (References 7, 8, and 9). 

At this point in the analysis, Equations 2-1, 2-2, and 2-4 represent three different formula- 
tions of a single problem, and each is a second-order differential equation with non-constant coef- 
ficients. The radius vector and true anomaly of a particle in an elliptic orbit in an inverse- square 
force field can be expressed as a function of the mean anomaly. Moulton (Reference 19) gives these 
functions as 


r 


a 1 - 


e cos M - ~ 2 ~ ( cos 2 M - 1 ) 


and 


v - M + 2e sinM + 4 e 2 sin 2 M , 


and, for small e, 


fji e r 3 - M e 3 jjL ~ e cos M - -g- (cos 2M - 1 ) - * * * 

- M e a ~ 3 [1 + de cos M + * * *] 

- 1 + 3e cos t 


and 


v - - 2e sin t , 


since M = nt andn = a" 3 = 1. 
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Thus, to first order in e, Equations 2-1 and 2-2 become 


d 2 <f> , N 

+ 3 (r 2 - rj (1 +3e cost)0 


2e sin t 


(2-5) 


and 


d 2 <£' / X 

-^ 2 - + 3(r 2 - r x ) (1 +3e cos t) <£ 


6e(r 2 - r 1 ) sin t 


Similarly, to first order in e, 


3 { r 2 " r ! ) + e COS V r -1 , 

1+ecosv = [3(r 2 - Tj) +e cos vj [l- ecos V + e 2 cos 2 V- ■ • -J 

- 3 ( r 2 " r ,) + [l " 3 (r 2 - rj )] ecos 


V , 


(2-6) 


and Equation 2-4 becomes 


d 2 <1> 

dT 2 



rj + [l-3(r 2 



e cos v f <$> 


2e sin v . 


(2-7) 


Reduction to Canonical Mathieu Form 

Each of the three formulations of the preceding section leads to a nonhomogeneous, linear, dif- 
ferential equation with periodic coefficients. Although in general there is no analytic method for 
obtaining solutions to this broad class of equations, a number of particular cases have been ex- 
haustively studied. One of the better known examples is Mathieu's equation. It is possible to trans- 
form each of Equations 2-5, 2-6, and 2-7 into the canonical form of Mathieu's equation by a suitable 
change of variables. The solution to this equation is well known and gives the solution to the pitch- 
motion problem. 

Historically, the equation was first studied by Mathieu (Reference 20) in conjunction with the 
vibration of elliptical membranes. McLachlan (Reference 21) has reviewed the general theory of 
Mathieu's equation, and its application to the analysis of loud speakers, frequency modulation, prop- 
agation in wave guides, oscillatory systems with periodic disturbances, etc. A point of interest is 
that the linear variational equation used in evaluating the infinitesimal stability of solutions to non- 
linear differential equations frequently has the form of Mathieu's equation; see Stoker (Refer- 
ence 22) and Cesari (Reference 14). 
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McLachlan gives the canonical form of Mathieu's equation as 


d 2 y 

+ (a™ 2q cos 2z)y - 0 . 


( 2 - 8 ) 


The substitutions 


2z = t , 
a = 12(r a - ri ) , 


and 


q = - 18 (r 2 - r t ) e 


transform Equations 2-5 and 2-6 to 


d 2 4> 

“7^7 + (a~ 2q cos 2z)0 - 8e sin2z 

dz z 


(2-9) 


and 


d 2 

-r + (a~ 2qcos 2z)q^' 

dz^ 


24e (r 2 - r 1 } sin 2z 


(2-10) 


The substitutions 


2z v , 

a = 12(r 2 -r,) , 


and 


2q = ( a - 4) e 


transform Equation 2-7 to 


d 2 4> 

- — r + (a - 2q cos 2z ) 4> - 8e sin 2z . 

dz * 


( 2 - 11 ) 
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The effects of the choices of independent variable and reference system could be determined 
even before these transformations were completed. A comparison of Equations 2-5 and 2-7 in- 
dicates that the choice of independent variable affects the periodic part of the homogeneous equa- 
tion but not the forcing function. The homogeneous equations differ for two reasons: first, the de- 
pendent variable is expressed as a function of v instead of t, and second, the dependent variable is 
changed from to a> by a time-varying transformation. Similarly, comparison of Equations 2-5 and 
2-6 indicates that the choice of reference frames affects the forcing term but not the homogeneous 
equation. In Equation 2-5 the forcing term is a result of the acceleration of the reference system, 
while in 2-6 it is due to the torque which would exist on a perfectly oriented (<£' = 0) spacecraft 
since the mean vertical is not a zero-torque axis. 

When the eccentricity, e, is zero, the true and mean anomalies are identical and there is no 
difference in the formulas— Equations 2-9, 2-10, and 2-11— resulting from the three approaches. 
There is an interesting difference when r 2 ~r 1 = 1/3, i.e., a = 4, for any non-zero value of e. The 
homogeneous part of Equations 2-9 and 2-10 have an unstable solution in the neighborhood of the 
line a = 4, as shown in the next section. However, Equation 2-11 is reduced to a constant coeffi- 
cient harmonic oscillator when a = 4, since q = 0 in that case. The difference is caused not by 
the change of independent variables from t to v, but by the change of dependent variable from <f> to 
0. The latter transformation removes, to the first order in e, the resonant term corresponding to 
1 + e cos v that causes Equations 2-9 and 2-10 to be unstable. 

Stability 

The complete solution of a general second-order, linear, nonhomogeneous differential equation 
can be written as the sum of two linearly independent solutions of the homogeneous form of the 
equation and any particular solution of the nonhomogeneous equation. The homogeneous part of 
Equations 2-9, 2-10, and 2-11 is identical with the canonical form of Mathieu T s equations; thus the 
nature of the complementary solutions is well known. 

There is a general result from Floquet theory (Struble, Reference 13) which says that any solu- 
tion, y, of a linear differential equation with continuous periodic coefficients of period r can be 
represented as 


y(z) ~ P(z)exp(/xz) 

where P(z) is periodic with period r. If y(z) is a solution of Equation 2-8, so is y(-z); thus the 
complete solution can be expressed as 

y(zj = Cj P(z)exp (/zz) + c 2 P(-z ) exp (-/zz ) . (2-12) 

It is apparent from Equation 2-12 that the solution to Equation 2-8 is stable if and only if ^ is 
a pure imaginary number. The form of /z, real, imaginary, or complex, is a function of the parameters 
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a and q. It is possible to plot the regions in which m is real or imaginary as a function of these 
parameters; such a plot indicates regions in which the solutions to Equation 2-8 are bounded and 
those in which the solution is unbounded. 

The stability of the complete solution of any of Equations 2-9, 2-10, or 2-11 requires that both 
the complementary solution and the particular solution of the equation are bounded. That is, for 
stable motion it is necessary but not sufficient that the parameters a and q yield a bounded solu- 
tion to Equation 2-8. 

The range of the parameters a and q for this application is restricted by the physical nature 
of the problem. Only spacecraft whose inertia range is in the Lagrange region will be considered; 
thus, 0 < r 1 < r 2 < 1. The orbital eccentricity must lie in the range 0 < e < 1 for a closed orbit; how- 
ever, it would be wrong to use these bounds in determining the range for q. All three formulations 
contain series expansions in which terms of second order and higher in e are dropped; this implies 
e « 1. If e « 1 is taken as e < 0.1, the range of q can be limited to -1.8 < q < 0 for the first two for- 
mulations and -0.2 < q < 0.4 for the last one. A second limit exists, since when e exceeds 0.6627 • • • 
the series developments for r and v used in formulating Equations 2-9 and 2-10 are not convergent, 
although the expansion for Equation 2-11 is not so restricted. Thus a range of 0 ^ a ^12 for all 
three equations, -12 < q < 0 for Equations 2-9 and 2-10, and -2 < q < 4 for Equation 2-11 completely 
covers the admissible range of parameters. 

Figure 2-3 shows the regions of bounded and unbounded solutions, mapped from the numerical 
data in Reference 23.* The regions are symmetric about the a axis; thus it is sufficient to plot 
only the absolute value of q to represent the range determined above. In using this figure to eval- 
uate the stability of the solutions of Equations 2-9, 2-10, or 2-11 it is necessary to convert the 
parameters r 2 - r x and e to the parameters a and q by the appropriate transformation. The con- 
version from ( r 2 - r 1 , e) to ( q , a) is not the same for Equation 2-11 as it is for Equations 2-9 and 
2-10, and neither of the sets of conversions is linear in the variables r 2 - r 1 and e. These two facts 
make it difficult to "see” the effects of changes in the physical parameters r 2 - r x and e, and obscure 
the differences resulting from the choice of independent variables. To avoid these difficulties the 
regions of bounded and unbounded motion are also demarcated with respect to the physical param- 
eters r 2 - r 1 and e in Figures 2-4 and 2-5. Figure 2-4 is applicable to both formulations that re- 
tain time as the independent variable, (i.e., Equations 2-9 and 2-10); Figure 2-5 corresponds to the 
choice of true anomaly for independent variable (i.e., Equation 2-11). 

These figures differ in two respects: first, there is a region of unbounded motion beginning 
at r 2 - r j = 1/3, e = 0 in Figure 2-4 but not in Figure 2-5; second, the boundaries of the unstable 
regions beginning at r 2 -r 1 = 1/12, e = 0 and r 2 - rq =3/4, e = 0 differ markedly between the 
two figures as e increases. In the previous section, the first of the two differences was attributed 
to the transformation <i>(v) = <£(v) (l + e cos v) which removes, to the first order in e, the resonance 
term that produces this region of instability. The instability is, of course, still present, but it can- 
not be shown with a solution to the first order in e such as is determined herein. The second 


*The lines on the figure should be thought of as lines of demarcation, for regions where a given value (q, a) is unstable. 
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Figure 2-3— Regions of instability for Mathieu's equation. 
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Figure 2-5— Regions of instability for Equation 2-1 1 . 
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difference is a more obvious result of discarding higher-order terms. When the boundaries are 
compared for e « 1 they are quite similar. 

It has already been noted that it is necessary to consider also the particular solutions of Equa- 
tions 2-9 through 2-11 before deciding that a given combination of parameters is stable. The im- 
portance of this can be seen by considering the case r 2 - r 1 = 1/3, 0 < e < 1, which is frequently 
called the "linear pitch resonance case." With these parameters. Equation 2-11 becomes 

d 2 0 

— — y +40 ~ 8e sin2z , 

dz z 


which is the equation of a linear oscillator forced at its resonant frequency. Thus the solution of 
Equation 2-11 is clearly unstable when r 2 -r 1 = 1/3 for any non-zero value of eccentricity. This 
is indicated by a dashed line in Figure 2-5. Particular solutions of Equations 2-9 and 2-10 may 
also be unstable in this region but, since the complementary solution is already known to be un- 
bounded, it is not necessary to pursue this aspect of the problem. 

Figures 2-4 and 2-5 and the above discussion show that there are three regions of unstable 
pitch oscillation for any non-zero eccentricity, however small. These regions reduce to values of 
r 2 - i-j equal to 1/12, 1/3, and 3/4 as e approaches zero. When these values of r 2 - r l are sub- 
stituted in the linearized equation for pitch motion in a circular orbit, it can be seen that r - r x 
= 1/12 corresponds to a pitch libration with a natural frequency of one-half cycle per orbit. Sim- 
ilarly, r 2 ~ r i - 1/3 anc * r 2 " r i = 3/4 correspond to natural frequencies of one and three-halves 
cycles per orbit, respectively. 

Damping 

The question of stability can be determined for a rigid spacecraft by knowing whether the value 
of f i corresponding to a given set of parameters (q, a) is real or imaginary. A specific set of 
parameters actually determines the numerical value of n and thus the rate of divergence of the 
pitch motion, when the point (q, a) is in an unstable zone. This rate of divergence becomes im- 
portant when the effects of energy dissipation (i.e., damping) are considered. 

A rigid spacecraft moving in an inverse- square gravitational force field is not subject to damp- 
ing. In order to model damping correctly it would be necessary either to consider a multiple-body 
spacecraft or to include other sources of torque such as reaction with the magnetic field. This is 
beyond the scope of this dissertation; however an idealized velocity- dependent damping term can be 
inserted in Equation 2-1, without physical justification and solely to obtain a "feeling" for the effect 
of damping. Thus, we consider the equation 


d 2 4 > d<£ 

IT* + 2K 3F + 7T 


( r 2 - r i) ^ ' 


(2-13) 
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When the expansions and simplifications of the second and third sections of this report are applied 
to Equation 2-13, the final result (which is analogous to Equation 2-9) becomes 


d 2 4> d<p 

~^~2 + 4 K dz + ( a “ 2q cos 2 z)<p - 8e sin 2z , 


(2-14) 


where, as before, 


0 ~ pitch angle , 

z = 2t , 
a = 12(r a - ri ) , 


and 


q - " 18 (r 2 -rje . 

The first derivative term can be removed from Equation 2-14 with the substitution 


y(z) - 0(z) exp (2k z) 


(2-15) 


which converts Equation 2-14 to 


d 2 y 
dz 2 


( a - 4k 2 - 2q cos 2z) y 


8e sin 2z 


(2-16) 


In this case the stability of the complementary solution to Equation 2-16 can be ascertained from 
the nature of the region surrounding the point (q, a - 4k 2 ) in Figure 2-3. Once the solution y(z) is 
known, the pitch motion is determined from the inverse of Equation 2-15 as 

0(z) = y(z) exp (~2kz) . 

Thus, if y(z) is in a stable region, then the complementary pitch solution is asymptotically stable. 
However, when the complementary solution to Equation 2-16 is divergent, it does not necessarily 
follow that the pitch motion is also divergent. In this case the solution to Equation 2-15 has the form 

0(z) = exp (-2 kz) [c t P(z) exp (/xz) + Cj P(z) exp (-/xz)j 

= c : P(z) exp [(/x- 2k)z] + c 2 P(-z) exp [-(/x + 2 k)z] 
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The pitch motion is therefore asymptotically stable, with damping, whenever Ul < 2 k . Values of ^ 
are mapped with regard to e and [3 (r 2 - r x ) - * 2 J 1/2 in Figures 2-6, 2-7, and 2-8, which are based 
on Reference 24. These charts differ from Figure 2-4 in two respects; first, each chart is ap- 
plicable only to a portion of the parameter space; second, the vertical coordinate represents the 
frequency of oscillation of the solution to Equation 2-13 for a circular orbit. 

The process of determining stability for the complementary solution to Equation 2-14 can now 
be reduced to an automatic process. First, given the parameters r 19 r 2 , e, and k 9 find j^3(r 2 -rj -* 2 ] 172 
then, determine the value of ^ from the appropriate figure for the point (e, [3(r 2 - r x ) - * 2 J 1/2 ^. The 
complementary solution is asymptotically stable if ^ < 2 x, periodic if p = 2 *, divergent if i± > 2 k. 

Thus the complementary solution to the in-plane pitch-motion equation can be asymptotically stable 
even for those combinations of e and r 2 - r x that are divergent according to Figure 2-4. Figures 2-6 
through 2-8 show that p increases directly with increasing e, and thus the amount of damping needed 
for stable motion also increases with increasing orbital eccentricity. 

The complete solution to Equation 2-14 consists of the sum of the complementary solution and 
any particular solution of the nonhomo geneous equation. Although the above analysis indicates that 
all of the regions of instability previously found for this equation can be eliminated by damping, it 
does not guarantee that the particular solution introduces no instabilities. However, the particular 
solution to this equation is analogous to the forced motion of a damped harmonic oscillator, and, as 
is shown in the next chapter, the forced response (or particular solution) is a bounded, slightly 
perturbed, sinusoidal oscillation, except when r 2 - r x = 1/3. 

Conclusions 

The motion of a gravity -gradient- stabilized spacecraft in an elliptic orbit differs in kind from 
the motion of the same spacecraft in a circular orbit. It is not correct to assume that the only im- 
portant effect of eccentricity is to add a forcing term to the constant- coefficient equations for a 
circular orbit. Actually, it is also necessary to include the periodic variation in the gravity- 
gradient restoring torque caused by the change in orbit radius with position for an elliptic orbit. 

A rigid gravity- gradient spacecraft moving in a circular orbit has a stable position of equilib- 
rium, when the principal axis of least inertia is aligned with the local vertical, and the principal 
axis of maximum inertia is aligned with the orbit normal. This can be expressed in the notation of 
this dissertation by a requirement that 0 < r x < r 2 < 1; the portion of parameter space satisfying 
this requirement is called the Lagrange region. When an elliptic orbit is considered, this stability 
criterion is no longer sufficient. In fact, three areas of parameter space in the Lagrange region 
are unstable in pitch even for very small eccentricities. As eccentricity approaches zero, these 
three regions correspond to a spacecraft with a libration frequency for a circular orbit of one- 
half, one, or three-halves cycles per orbit, respectively. 

The equation describing in-plane pitch motion can be developed in a number of ways. The 
formulations of Beletskii (Reference 8) and Baker (Reference 4) are compared with a third formu- 
lation that has time for an independent variable, as did Baker, but relates the motion to the local 
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vertical, as did Beletskii. The regions of unbounded pitch motion are mapped with regard to the 
inertia parameter, r 2 - r x , and orbital eccentricity, e, for each of the formulations. All three ap- 
proaches lead to equations that are unstable when r 2 -r 1 = 1/12, 1/3, and 3/4, for all non-zero 
values of e. The unstable regions expand in size with increasing eccentricity; Figures 2-4 and 2-5 
map the regions of unstable motion with regard to eccentricity, e, and the inertia difference, r 2 - r 

When the spacecraft parameters and orbital eccentricity are in one of the unstable regions, the 
pitch motion consists of exponentially increasing quasi-periodic oscillations. The rate of expo- 
nential growth, fx 9 is determined by the parameters r 2 - r x and e; Figures 2-6 through 2-8 plot lines 
of constant ^ with regard to these parameters. Energy dissipation in the form of velocity-dependent 
damping of sufficient magnitude can suffice to eliminate all the instabilities discussed above. The 
magnitude of the damping coefficient, *, required to stabilize the motion of a particular spacecraft 
increases directly with increasing values of e. 




CHAPTER 3 


PLANAR PITCH LIBRATION-ASYMPTOTIC EXPANSION THEORY 


Introduction 

There are several approaches to analyzing a motion when the equations describing it are not 
easily solved. The most common approach is to consider an equation which is T, nearly the same 
as” the equation of interest but whose solution is known. The suitability of this approach is de- 
termined by the nature of the approximations required to reduce the original equation to the form 
of an equation with a known solution. In Chapter 2 it was necessary to discard terms of order e 2 
and higher in the expansions of r and v as functions of t, and for the most part only the solution to 
the homogeneous portion of the equation was discussed. 

A second approach is to obtain an approximate solution to the specific equation being studied. 
This chapter presents a solution to the in-plane pitch -libration problem in terms of asymptotic 
expansions in powers of a small parameter (e). This method permits inclusion of terms to any 
desired order in e (although it will be shown that inclusion of terms of order e 4 and higher does 
not introduce new phenomena), and leads to the development of a particular solution of the non- 
homogeneous equation as well as inequalities defining regions where the solution of the homogene- 
ous equation is unbounded. 

An asymptotic solution to the in-plane pitch problem with true anomaly as the independent 
variable was developed through the second approximation by Beletskii (Reference 8). He obtained 
a particular solution of the nonhomogeneous equation, similar to the first terms of the solution 
found herein (subject to the difference in variables). He also noted the existence of the first 
resonant zone, and the linear resonance for the particular solution derived herein. However, in 
his formulation it is necessary to consider the third -order solution (which has not been done) to 
find the second and third resonant zones developed in this chapter. 

Equation of Motion 

This chapter, also, is concerned with the linearized equation representing in-plane pitch mo- 
tion of a rigid gravity -gradient spacecraft moving in an elliptic orbit. The equation is 


d 2 4 > 
dt 2 


3 ( r 2 ~ r l) 71 * 


+ 


29 


d 2 V 
dt 2 


(3-1) 



where 


and 


r = r(t ) 


v - v(t) 


Expansions of (a/ r) 3 and v-M in series involving powers of e and the mean anomaly are given by 
Cayley (Reference 25) through terms of order e 7 . The first terms in these expansions are: 


(f 


3 e 2 15 e 4 


1 + 9 + 


9 e ' 


g + 3 e cos t + ~Y~ cos 2 t 


‘27 

53 

4 [7 77 1 

g cos t 

+ "g“ cos 3 t 

+ 2 cos 2 t + g cos 4 t| 


and 


v - M + 2 e sin t + 4 e 2 sin 2 t + e : 


r i 13 

- 4 sin t + ^2 sin 3 t 


+ 


f 11 103 ] 

- 24 sin 2t + ~g^~ sm 4 t| 


With the aid of these two expansions, and with the substitutions 


3 ( r 2- r i) 


U 3 2 = 3 ( r 2' r i) 


3 e 2 15 e 4 


1 + — + 8 


Equation 3-1 can be written as 

d 2 <p ! „ ^ „ I_ 9e 2 


+ <gj 2 +co 2 
dt 2 ] 3 3 


r 27 53 


3 e cos t + — 2~ cos 2 t + e 3 l“g- cos t + cos 3 t 


A (7 77 \ 

+ e 4 \~2 cos 2 t + ~g" cos 4 tJ 


4 > ~ 2esint + 5 e 2 sin 2 t 


e 3 . 39e 3 . _ lie 4 . _ 103e 4 . . , n nS 

sin t + — 4 — sm 3t - — g — sin 2t + — g — sin 4t . \o-&) 
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As in Chapter 2, this is a nonhomogen eous linear second-order differential equation with periodic 
coefficients. Its total solution is composed of any particular solution to the nonhomogen eous equa- 
tion plus the complementary solution to the associated homogeneous equation. 

Asymptotic Solutions 

Although Poincare wrote the initial discussion of asymptotic expansions in powers of a small 
parameter, the method used herein was developed by Krylov and Bogoliubov (Reference 26) and in 
more detail by Bogoliubov and Mitropolsky (Reference 27). Appendix C presents a particular set 
of formulas that specialize and extend the approach of Reference 28 to include third- and fourth- 
order terms for an equation of the form 

d ^ 

+ w 3 2 0 - e f j (0, t) + e 2 f 2 (0, t) + e 3 f 3 (0, t) + e 4 f 4 t) . (3-3) 

The development is routine; readers interested in the underlying theory or a more detailed dis- 
cussion should refer to Bogoliubov and Mitropolsky. 

In essence the method reduces the solution of Equation 3-3 to a set of simpler problems. Spe- 
cifically it seeks to find functions 

UiCaJ, t), A. (a) , B i (a) , i - 1, ■ * * , n 

such that 

0 - a cos v y + eu 1 (a,'P, t)+***+e m u m (a, 'P, t), (3 -4) 


where 


da 

dt 


e Aj (a) + 


+ e m A (a) 

m v ' 


and 


d<P 

dtT " ^ + e Bj (a) + 


+ e m ( a) 

m v ' 


satisfy Equation 3-3 to order e m+1 . When there is a resonance between the perturbing expressions 
and the natural frequency gj of the unperturbed system, the problem is restated to find functions 

u. (a, 0, t) , A. (a, 0) , B. (a, 0) , i = 1, * * * , m , 
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such that 


4> - a cos + ©) + e Uj (a, 0, t ) + * * * + e m u m (a, ®, t ) 


(3-5) 


where 


da 

dt 


e A l (a, 0) + * * • + e m A m (a, 0) 


and 


d© 

dt = e Bj (a. ®) + *** + e m B m (a,0) 

satisfy Equation 3-2 to order e m + 1 for parameters in which a 2 ^(p/q) 2 . These two formulations 
are referred to as the ’’nonresonance case” and the ’’resonance case,” respectively. 

The following section develops a solution to Equation 3-2 for the nonresonant case; it will 
serve as a particular solution to the nonhomogeneous equation. Examination of this solution will 
also reveal the values of a- at which resonance occurs. This is followed by a section that develops 
inequalities that determine when the resonant solution to the homogeneous equation associated 
with Equation 3-2 is unbounded. These inequalities are plotted in a form similar to that of Chap- 
ter 2. 

A Particular Solution to the Pitch-Libration Problem 

When Equation 3-2 is expressed in the form used in the method of asymptotic solutions, the 
result is 


d 2 


dt 


2 + oj£ <£> - e f j (<£, t) + e 2 f 2 (<£, t) + e 3 f 3 (<£, t) + e 4 f 4 (<£, t) 


where 


f t (</-, t ) 2 s in t - 3v 2 (cos t ) , 


f 2 (i,t) - 5 s in 2t - ~~ 2 ~ (cos 2t) 


f 3 (0> t) - ~ "4 sin t + 


39 . (27a 2 53^ 2 _ \ , 

-4- sin 3t - l — g — cos t + g cos 3t ) <f> 


(3-6) 
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and 


f 4 (<*. t) 


11 . 103 /% 2 „ 77Z 2 \ 

g- sin 2t + -g— s in 4t - 1 ~ y cos 2t + g~“ cos 4t\p. 


The formulas which are used in obtaining asymptotic solutions to equations of this type are 
developed in Appendix C. In essence Equation 3-4 is substituted in both sides of Equation 3-6 and 
terms of the same magnitude in e are collected. The first order terms give 


^5 2 u 1 (9 2 u 1 d 2 u t 

’ 2 + 'dWit + 77T + 


Fj + 2^A 1 sinW + 2acoB 1 cos V 


where 


F t - f ! (a cos 'P, t ) 

and as explained in the appendix this equation can be solved for A 1? B x , and u r These results can 
then be used to expand similar equations for f 2 , A 2 , b 2 and u 2 and the process can be continued 
to any desired order. The computations are rather tedious and only a few of the intermediate re- 
sults are included in the remainder of this section. 

In the first approximation, 


2 s in t - 


3au 


[cos ( t + V) + cos ( t - *P)J , 


A! (a) - 0, 

B x (a) - 0, 


and 


Uj (a, W, t ) 


2 s in t 3acc 2 cos ( t + M / ) cos ( t - W) 

+ 2- 1 + 27^- + r--^T_ ' 


The first -order solution to Equation 3-6 is 

<$> ~ a cos W , (3-7) 
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where 


da 

dt 


cM 

dt ’ ^3 * 


In the second approximation 


F 0 

r 2 


^2 - 

- 1 


sin 2t - 


9a£ 4 9aa: 2 

cos 4 / - 


2(l-4a 3 2 ) 


(i + rrt-) <2t + v> + (i + 


cos (2t - W) 


a 2 0 0 , 


9a 4 


4-3 (1 " **i) 


and 



L 1 

sin 2t 9acr 2 

1 . - 2 ^ 

cos ( 2t — ) /. £ 2 \ cos ( 2t - *P) 

U 2 = 

l 

CM _ 
PO 

ci 

1 

!*l 

-,*-4 ’ 16 | 

V 1 1 + 2-3 J 

1 + ^ 3 + l 1 + 1 - 2 cc 3 j 1 -a.-3 


The second-order solution to Equation 3-6 is thus 


£ 


a cos ¥ + 


3ae£ 2 

2 


cos ( t + 'P) cos ( t - ¥) 
1 + 2^3 + 1 - 2-3 



sin t 


(3-8) 


where 


da 

dt 


0 


and 


dJP _ 9e 2 £ 4 

dt = + 4oj 3 (1 - 4oj 3 2 ) ' 
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In the third approximation 


?3 g t sin t + g 2 sin 3t - a [g 3 cos (t + V) + g 4 cos (t -<P) + g 5 cos (3t +*P) + g 6 cos (3t -<P)] 


A, = 0, 


B 3 = 0, 


and 


gj sin t g 2 sin 3t 
“• = * * 


where 


l ~g 3 cos (t + 

L 1 + 2“3 


¥) g 4 cos (t -<P) g 5 cos (3t + ¥) g 6 cos (3t -<P)"| 


+ ■ 1 - 2^3 + 


. )' 

3" (3 + 2o) 3 ) + 3(3- 2o. 3 ) J - 


jl_ ^ _ & 2 r 35r 2 

4 o/, 2 _ i \ o/.. 2 a\ 5 2 1 

0: 3 1J 


2 (' 3 2 - 1 ) " % 3 2 - 4 ) 


, 39 _ 9£ 2 _ 3S: 2 r 3£ 2 1 

g 2 4 o I. . 2 i \ n / . 2 ,, \ 5 - ^ 2 _ jJ 


2(^3 *) 2 (^3 2 4) 


_ 27^ 27S 4 / S 2 \ 27S 4 

g 3 ‘ 16 + 32(1 +o; 3 ) l 1 + 1 + 2cc 3 J + 8 (1 - 2^ 


27S 6 


) " 4o.' 3 (l - 4 c<; 3 2 ) 


/"3 + 1 \ 

V 1+2 -3/ 


27S 2 

*4 “ 16 + 


27 S 4 / 2 r 2 \ 27£ 4 27<£ 6 (" 3 ~ 1 \ 

32 (l-. 3 ) + l-^ 3 j + 8(l + ^ 3 ) 4 x- 3 (l - 4 g, 3 2 ) 


53S 2 27cc 4 ( 2 2 \ 27£ 4 

16 + 32 (1 + o, 3 ) ^ + 1 t 2aJ + SflTi^J 


5 3£ 2 27S 4 ( S 2 \ _2_7£ 4 _ 

g 6 - 16 + 32 (1 - o; 3 ) ^ + 1 - 2w 3 y + 8 (1 - 2o; 3 ) ’ 

The third -order solution to Equation 3-6 becomes 


, _ 3aeS 2 

<P - a cos ^ + 2 


ii [cos (t +«P) cos(t-V)" 
1 + 2o; 3 + 1 - 2co 3 


9ae 2 co 2 
16 


, — y- 2 -Y cos (2t + ¥) ( 7? \ cos (21:-^) 2e e 2 / 3S 2 \ 

[_V 1 + *v i + «, + V + r^j-(i -^ + ^i sint + ^7(f-Z^rr> 
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where 


and 


da 

dt 


0 


dV 

dt 


9e 2 Z 4 


4 , 3 (1 - 4a- 3 2 ) 


As the order increases it becomes harder to find a solution. However, examination of each of 
the three preceeding solutions —Equations 3-7, 3-8, and 3-9— indicates that each new approximation 
introduces a significant new result. In particular, u x has an unbounded term when = 1/2, u 2 
when < 3 = 1, and u 3 when = 3/2. However, further approximations do not produce further 
meaningful discrepancies. In particular, the fourth approximation introduces new terms of the 
form 

h x cos (4t +r V) h 2 cos (4t -*P) 
u 4 T ae “ 4 8(2T7,^) + 8(2-T 3 j~ “ 

+ fourth-order corrections to terms found before, 


2 - 16 


s in 4t 


and 


A 4 - 0 , 

B 4 * e 4 h 4 , 

where h. (i = 1, • • • , 4) are constants that could be evaluated with Appendix B formulas. How- 
ever, since 0 < ^ 3 2 <3 for small e, the terms in u 4 are bounded and there is no need to consider 
them further. 

Instead, return to the third-order approximations. Since da'dt = 0 for all a, a particular solu- 
tion to Equation 3-6 is 


<t> = 


2e 

T7i Slnt 

i 1 




Si 


sin t 


4 



sin 3t 


(3-10) 


This solution represents the motion of a spacecraft for the particular case in which a 0 = 0. It 
also approximately represents the forced motion of a rigid spacecraft in an eccentric orbit when 
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"light damping" is present. Equation 3-10 is a particular solution to the nonhomogeneous Equa- 
tion 3-6. In the next section the solution to the homogeneous equation associated with Equation 3-6 
will be found for each of the regions 1/2 , 6j 3 ^1 ? and oj 3 ^ 3/2. The combination of these com- 
plementary solutions and the particular solution, Equation 3-10 is sufficient to determine completely 
the nature of the solution to Equation 3-6 in the region of interest. 


Regions of Resonant Pitch Motion 

An asymptotic solution to the in -plane pitch-libration problem was developed for the non- 
resonance case in the previous section. This solution represented unstable motion for a> 3 = 1/2, 

1, or 3/2, which indicates that a resonance-case type of solution such as Equation 3-5 should be 
used in these regions. The process is quite similar to that of the preceding section, in that solu- 
tions are built up with terms of increasing order of approximation. 

It was shown in the previous section that the first three approximations introduce all the reson- 
ant phenomena of interest for this problem. Fourth-order terms are therefore ignored in this 
dissertation. 

The homogeneous part of Equation 3-2 can be written to the third order in e as 

d 2 l 

+ a, 2 <f> - efj (/., t) + e 2 f 2 O, t ) + e 3 f 3 (*, t ) , (3-11) 


where 


fj - - 3^ 2 (cos t)0 , 


9£ 2 

f 2 ~ “2 ( cos 2t ) 4 , 


and 


f 


3 


21a: 2 

8 


(cos t ) {f. 


5&: 2 
' 8 


(cos 3t ) «£ . 


Solutions of Equation 3-11 will be developed for co 3 ^1/2, 1, and 3/2; as before, most of the algebra 
will be omitted. The first series of solutions will assume oj 3 ^ 1/2, and define ^ 3 2 - 1/4 = eA. In the 
first approximation, 


Fi 


- 3aai 2 

2 


[cos ( t + l P) + cos ( t ~ ^)] , 
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3oj 2 

2 a sin 20 , 


3£ 2 

2 


cos 20 + A , 


and 


u 


i 


3a£ 2 

4 


cos ( t + V P) . 


Then, in the first approximation, 


where 



(3-12) 


da _ 3e^ 2 

2 a s ^ n 2© 


(3-13) 


d© 

dt 


3ea 2 

2 


cos 20 + 


4 ’ 


(3-14) 


and a is not a constant. When new variables u = a cos©, v = a sin© are introduced into Equation 
3-13 and 3-14, the system 



results. This can be reduced to either 


or 



o , 



0 . 
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Both u and v are represented by unbounded oscillations when 


i , ll 

3e£ 2 

K 4 1 < 

to 


(3-15) 


and, since a ~ [u 2 + v 2 ] 1/2 , inequality 3-15 is also a condition for unbounded oscillation of <£. 

A more exact boundary can be found by considering the second approximation, in which 


9a£ 4 3aA^ 2 , . 9acu 2 

g cos V + 2 cos ( t + *P) - 4 


3^2 

1 + ~ 2 ) COS (2t + V) + II - —) cos (2t + *P) 


A = 0 , 


-A 2 + 


27^ 


and 


■t a a cos 


(2 + ®) 


3eaa 2 


, 3t 

COS [“o’ + © 


(3-16) 


where 


da 3eoj 2 

dt = 2 a sin 20 


d© 


3ecc 


dt - ^2~ cos 20 + K 2 - 4 


-i - 4 


l\ 2 27e 2 


8 


The same process that was used in the first approximation indicates that Equation 3-16 represents 
unbounded motion whenever 


4 - ("3 2 ' 4 ) 


\ 2 27e 2 £ 4 

3eS 2 

■+ 

00 

< I 2 


(3-17) 


It is possible to refine the boundary of the unstable region near a> 3 = 1/2 still further by consider- 
ing higher -order approximations, but the gain is small. 

The unbounded motion for eo 3 * 1 appeared in the second -order approximation (Equation 3-8) to 
Equation 3-2 and will likewise only be apparent in the second-order resonant solution. In this case 
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co 3 is nearly one; i.e., eA s a 3 2 - 1, and, in the first approximation. 


„ 3a2; 2 „ 

- - 2 [cos (t + *P) + cos (t - ¥)] . 

A 1 - 0 , 


B 


l 


A 

2 • 


u 


l 



COS ( t + V) - 


3ao: 2 

2 cos ( t - V) , 


and 


<t* ~ a cos ( t + ©) , 


where 


da 

dt = °- 


d© (c i * 
~dt 2 


Thus, in the first approximation the amplitude of the oscillations remains constant. This is not the 
case when the second approximation is considered, for 


F 


2 


4 aAy 2 3a\7 2 

2 cos T + 2 cos (t + V) + 2 cos (t - V) 


9aa; 2 

4 



cos (2t + *F) - 


9a£ 2 

4 


(l - 6^ 2 ) cos (2t - *P) , 


9<£ 


g (l - a; 2 ) a sin 20 


(»-; 


} cos 2© - g 


3A 4 

4 


a AS 2 

- g cos (t + W) + 


3aAx 2 

2 


cos ( t - ¥) + 


9a£ 2 

32 



cos (2t + V) , 
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and 


^'2 

aeoi* 

0 - a cos ( t + 0 ) + 2 cos ( 2t + 0 ) . 


(3-18) 


where 


da 9e 2 d ; 2 ^ 

dt ~ 8 (1 - oj 2 ) a s in 20 , 


d 0 

dt 


«3 2 i 1 K-jT 


3e 2 d; 4 9e 2 S 2 , ^ v 

4 + g (l - oj 2 ) cos 20 . 


This set of equations represents unbounded motion if 


K~l ) 2 _ 3 e = 


(1- - 2 )l 


(3-19) 


This boundary can also be more precisely defined by considering the third approximation, in which 


F 3 - 2aAd '- 4 cos X V 


27ad 2 


16 


~27slo. 2 


9ad : 4 ( N 

7aA 2 

£ 2 3a^ 6 ~| 


16 

(1 + 2 a 2 ) + 

64 ( 3 + ^ ) 

+ '24 

+ 2 I cos 

(t +¥) 

9 ad 4 

15aA 2 3 2 " 

1 cos ( t - M') + 

a/Vl 4 

3aAd- 2 . ."1 


+ 8 

f ' 8 



16 (3 + £ 2 ) 

cos ( 2 t +M») 


9 a Ad 4 

4 cos ( 2t - X V) , 


9/Vd 4 

A 3 g a sin 20 , 


A 3 ^ 11A2'4 9Ad^*_ 

16 + 8 + 8 cos 20 , 


and 


0 - a cos ( t + 0 ) 


3aeA 2 . 9 > aed; 2 , 

+ 2 -(ij J - 1) COS0+ ■ -g- (4-< 


2 ) cos ( 2 t + 0 ) 


9 ae^ 

32 


2 £ 2 ( sA 

12 ^ + 3 J 


cos (3t + 0 ), (3-20) 
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where 


da 9e 2 qj 2 

dF " 8 


[1 -s 2 +s 2 


K 2 - !)] 


a sin 2© , 


d© 9e 2 S 2 
dt " 8 


'K 2 - 1 )] 


cos 20 + 


1 ~ *) 2 K 2 ~ O 3 3e 2 S> 4 lie 2 S 4 


K a - 0 


The motion represented by Equation 3-20 in the regions *1 is unbounded when 


K-l) 2 (-3 2 - 1 ) 3 3e 2 


2 4 1 le 2 S 4 


( a 3 " X ) 


[l -£ 2 + 2 ! 2 (cc ' 3 2 - i)] 


(3-21) 


The last region to be considered is in the vicinity ^ 3 %3/2. The existence of this unbounded 
motion was indicated only in the third-order approximation solution, Equation 3-9, and that is also 
true when the resonant solution is developed. In this case eA = cu£ - 9/4 and, in the first 
approximation, 


3aa> 2 

2 [ cos ( t + V) + cos ( t - ¥)] 


Aj = 0, 


- 3 


3aC; 2 3a oj 2 

Uj g cos(t+ l P)- cos(t- l P) 


- a cos 


(t * 9 ) 


where 
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In the second approximation, 


_ 9 a £ 4 a /£ 2 a /£ 2 

F 2 - 16 cos ¥ + — 4 — COS (t + MO + — Y~ cos (t “ V) 


9a£ 4 

4 



cos ( 2 t +W) 


9a£ 4 

4 



cos ( 2 1 + V) 


a 2 = 0 , 


B . Al 3^ 

B 2 " 27 16 1 


a/£^_ a/VT 2 9a .^ 2 ( Z 2 \ 

16 cos (t + ^) + 4 cos (t - V) + ~ 4 o~ [l + ~d7J cos (2t + V) 


9 ag' 2 

8 



cos ( 2 t - MO . 


and 



(3t \ 

3ae £ 2 

0 ( t 

/St 

\" 

<£ = 

a cos + 01 

~ 8 

2 cos ( 2 + 0 

“ cos ^ 2 + ( 

v\ 


where 


and 


da 

dt 


0 


d0 ^3 2 “4 H*" - ?) 3e 2 £ 4 

dt " 3 27 " ‘ 16 * 
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Thus even in the second approximation there is no sign of the resonance. However, in the third 
approximation, 
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where 


da 

53e 3 

CO 2 

9e 3 

L 


a sin 20 , 

dt “ 

48 


8 

l 1 

“ 4/ 


L_ 



\ 



d© 

53e 3 

& 

, 9e 3 oi 4 

L 


“3 2 -i ( w 3 2 -t) 4a 2 - 1) 

dt " 

'48 


8 

l 1 

■ 4 ! 

cos 20 + 3 - 27 + 243 


3e 2 £ 4 1 3e 2 _Sl / 2 _9\ 

" " 16 + 96 \^3 4] 

This system of equations leads to unbounded oscillations whenever 


2 ± 

“ 4 

( 3 -V 

\"3 “ 4/ 

4’a 2 - ?) 3 

3 “ 

27 + 

243 


3e 2 ^ 13e 2 ^ 4 / 9_\ 

16 + 96 \^3 ~ 4/ 


53e 3 co 2 

9e 3 £ 4 / 


" 48' " 

8 l 

1 ' 4/ 


(3-23) 


It would be necessary to consider the fourth-order approximation if a more precise boundary of 
this region of unstable motion was desired. 

The above analysis has shown that there are three regions in which the solution to the homo- 
geneous equation related to the linear pitch libration problem represents unstable motion. Figure 
3-1 demarcates these regions in which this occurs, with respect to r 2 - r 1 and e. 

Both this and the previous chapter have been concerned with the linear, in-plane, pitch motion 
of a rigid gravity -gradient spacecraft in an eccentric orbit. The most important result of both 
chapters is the demonstration that stability criteria developed for circular orbits cannot be ex- 
trapolated to even slightly eccentric elliptic orbits. 

It is known that the linear equation describing the attitude motion of a gravity-gradient- 
stabilized spacecraft has a stable equilibrium in the circular -orbit case when the inertia param- 
eters are in the Lagrange region, (i. e., when 0 < r x < r 2 < 1). Actually this region produces motion 
asymptotically stable at the equilibrium when there is T ’pervasive damping,” no matter how small 
the damping is (Pringle, Reference 16). Neither of these conclusions is valid for a non-circular 
orbit, no matter how small the eccentricity. In fact, there may be divergent oscillations for some 
inertia combinations even with light but pervasive damping. 


45 




These instabilities can be attributed to two causes. First, there is the forced response due to 
the non-uniform orbital rate. This enters the equations as a forcing term (Equation 3-1) and causes 
the equations of motion to be nonhomogeneous. Equation 3-10 shows the nature of the forced re- 
sponse; evidently the only important instability occurs when a > 3 = 1 (i.e., when r 2 - r t = 1/3). This 
same equation also indicates the steady-state behavior of a lightly damped spacecraft, except in the 
regions of parametric resonance. 

The second effect of orbital eccentricity is that the radius vector to the spacecraft, and thus 
the restoring torque, varies with time. This causes the equations of motion to have periodic coef- 
ficients and leads to a phenomenon called "parametric resonance." When cj 3 ~1/2, 1, or 3/2 (i.e., 
when r 2 - 1/12, 1/3, or 3/4), the amplitude of the pitch oscillation may increase without bound, 

to the order of approximation of the equations, even for arbitrarily small values of e. Figure 3-1 
is plotted for a series of inequalities that define regions of unbounded oscillation; see Equations 
3-17, 3-21, and 3-23. 

The preceeding two chapters have shown that it is not valid to extrapolate results from the 
linearized equations for a circular orbit to the linearized equations for an elliptic orbit. The re- 
mainder of the dissertation will demonstrate that there are inertia parameters in the Lagrange 
Region for which the linearized equations are not adequate for either circular or eccentric orbits. 
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CHAPTER 4 


GENERAL HAMILTONIAN EQUATIONS OF MOTION 


Introduction 

A complete description of the motion of a rigid gravity -gradient spacecraft in an eccentric 
orbit leads to a complex system of equations. Chapter 2 discusses some assumptions frequently 
used to simplify these equations. One of the most common is that it is valid to linearize the equa- 
tions and then treat the in -plane or pitch motion as independent of the coupled roll -yaw equations. 

A second assumption, shown to be false in Chapters 2 and 3, is that a slightly eccentric orbit does 
not affect stability except for the linear resonance case when the pitch period and orbital period 
are identical. In fact, the situation is more complex than the results of the previous chapters in- 
dicate, because linearization also leads to incorrect results even for small angle motion. 

The coupling between the pitch and roll-yaw equations has a significant effect on stability dur- 
ing large amplitude motion, as DeBra (Reference 6) showed by numerical integration. Kane (Ref- 
erence 11), using a Floquet analysis of an only partially linearized set of equations, showed that 
this coupling was important even for small -amplitude (e.g. 1 degree) motions. Breakwell and 
Pringle (Reference 12) confirmed this conclusion in a paper that treated the problem of nonlinear 
coupling with canonical transformations and the method of averaging. This very significant paper 
was the first to apply the methods of analytical mechanics to the three-degree-of -freedom problem, 
although Liu (Reference 28) had previously applied canonical transformations to the in-plane libra- 
tion of a spacecraft on an elliptic orbit. 

This chapter and Chapter 5 extend the method of Breakwell and Pringle to higher-order effects 
of orbital eccentricity on the motion of a rigid gravity -gradient -stabilized spacecraft in an eccentric 
orbit. The Hamiltonian for the motion is first developed to third order in terms of the pitch, roll, 
and yaw angles and their associated momentum variables. Then a series of canonical transforma- 
tions to new variables (for which the second -order portions of the Hamiltonian represent three 
uncoupled linear oscillators) is developed. The Hamilton- Jacobi equation is solved to obtain a 
third canonical transformation, which introduces cyclic coordinates in the second-order Hamiltonian. 
When the third rorder Hamiltonian is expressed in the new system, the resulting Hamiltonian equa- 
tions represent a set of perturbation equations that show the effects of nonlinear resonance. In 
Chapter 5 these perturbation equations are averaged and the inertia parameters that lead to various 
types of resonance are plotted. 
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Hamiltonian for Three Degree of Freedom Motion 

A development of the Hamiltonian for this problem follows in a straightforward manner from 
the kinetic and potential energies of the system. The kinetic energy, T, associated with the attitude 
motion of the spacecraft is 


T 


J. 

2 


( r i n i 2 


+ r 


:<V +n 3 2 ) 


where 

- \p ~ (c£+v)sin$, 

0 2 - c cos yp + (<£ + v) cos 0 sin p , 

fi r “ 0 sin 0 + (<p + v ) cos 0 cos yp . 

The potential energy V associated with attitude motion is 


V 



cos 2 <p cos 2 t' + (r 2 ~ l) (cos <p s in fc' cos \p + sin<£sini/-) 2 J 


and the Lagrangian L is 


L 3 T - V . 

The standard procedure at this point would be to introduce generalized momenta p., defined as 

<9L 


where 




^2 


and 


q 3 = <*■ • 
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Unfortunately this leads to complications when terms of similar orders of magnitude are collected. 
The general nature of the motion of a successfully stabilized spacecraft is well known. There will 
be small oscillations of 0 and 6 about zero, and a small oscillation of <j> about the local vertical 
which rotates at a rate of 2 tt radians/orbit. Thus the angular momentum about the ’Three” axis, 
p 3 , would have a steady component and a component representing small oscillations. However, 
this can be eliminated by defining 


dL 


dq 


3 


p 3 + I 3 n 


p 3 + 1 , 


with the happy result that all the q i and p. are small as compared with 1. The ’’order” of a given 
term should be understood as meaning the order of the products of the q.’s and p/s in the term, 
(e.g. q 2 2 p 3 is 3rd order in coordinates and momenta). 


dL 

90 


= rjflj. 


dL 

d'e 


- r Q cos 0 - fl sin0 


and 


p , + 1 = 


dL 

d<f> 


- r x sin 9 + r 2 ^2 cos & sin ^ + ^3 cos ^ cos ^ 


The Hamiltonian, H, is defined as 


H 


3 

dL . 

ak k qi 


- L(q, q, t) ; 


thus 


r i °1 2 +r 2 n 2 2 +n 3 2 i 


! J + j^r : 0 1 sin (9 - r 2 H 2 cos 6 sin 0 - 0 3 cos 9 cos 0 J v + V 


or 


H 


*1 


( r i n 0 


2 



( r * n 2 ) 2 



2 

3 


(p 3 + 1) v + V . 
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At this point the series expansions 


V 3 v 2 

s in 77 - 7) g- + *** , cos 77 - 1 - ~2 

are introduced into the Hamiltonian. When terms of equal order in the generalized coordinates 
and momenta are collected and the substitutions D q 2 D 0 and q 3 = <£ are introduced, the 
result becomes 


h = £> 


where 


v , 


Hi “ P3O-V) 


1 2 1 ,1,3m 

277 Pi + 2r~ P 2 + 2 ?3 + T 71 


(t) [( r 2 ~ r i h 3 2 + U ~ r i) h 2 2 ] 


1 il *2 

+ 2 r„ 


1 - r„ 


Q! 2 + 2 H 2 + r Hi P 2 + q 2 Pi 


1 - r 2 

H 3 = r, Hi P 3 (q t +P 2 ) + q 2 p 3 (q 2 + Pi) 


3_Vj_ 
2 a 3 


(t) t 1 ~ r 2 ) 


H 2 H 3 ’ 


and fourth-order and higher -order terms will be ignored. The series developments for r and v 
introduced in Chapter 3 will be used in the form 

v - 1 - T, (t) , 


where 


5e 


Tj (t) n ~ 2e cos t y~ cos 2t - [- cos t + 13 cos 3tj 


- ~ 24 ~ [“ 22 cos 2 t + 103 cos 4t] 


(4-1) 
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f 

1 

f 

i 

- and 


where 

(t) 3 = T + T a (t) 


, = - (. , 3e 2 15e 4 \ 

k 3 l 1 2 8/ 


and 


4 



T 2 = 

s i cos (it) 

(4-2) 



i = 1 


where 


27e 3 



S 1 

3e + 8 ' 




9e 2 7e 4 



S 2 

2 + 2 ’ 




53e 3 



^3 

8 > 


and 


77e 4 



S 4 “ 

8 


With these substitutions and with terms 

T 2 q 2 2 

and T 2 ^ 3 2 assigned to H 3 , H can be written 




3 



H 

kl 

1! 

(4-3) 



i = 0 


where 

Ho = 

T x (t) - 1 , 


\ 


53 


\ 

i 






Hj = P a T x (t) 


*2 " 2 


1 - r. 


— ' qi 2 + (l+k-kr 1 )q 2 2 + k(r 2 -r 1 )q 3 2 


+ 77 p i 2 + 77 p 2 2 + p 3 2 


1 r 2 

“77“ *1 P 2 + ^ 2 Pi 


and 


1 - r 2 3 

H 3 = ~7: ^1 P 3 (^1 +p 2 ) + ^2 P 3 (^2 + Pl) + 2 


k 

3 + T 2 ( t > 


) (l “ r 2 ) Qi 


q? q 3 


2 2 ) [( r 2 _ r 0 q 3 2 + i 1 ~ T 1 ) P 2 2 ] 


Diagonalization of Second Order Hamiltonian 

When Hamilton's equations are formed from Equation (4-3), the result is six equations that 
are nonlinear and coupled (i.e., terms q j3 p. j / i appear in the equations for q i or p.). The 
Hamilton equations have the form 


<9H 1 dH 2 (9H 3 

" dp { + dp { + ~dp i ’ 

and 

( ?H 1 dU 2 dH 3 

" Pi " dq. 1 Jq~ + <9q. 

Since n l contains only first-order terms, it cannot contribute any coupling terms to these two 
equations; however, both H 2 and H 3 may do so. If the H 3 term is ignored (i.e., if the equations are 
linearized) it is possible to find a canonical transformation to new coordinates for which the second- 
order Hamiltonian equations are uncoupled. It is then possible to use the solution to these equa- 
tions to obtain what amounts to a set of six perturbation equations that determine the effect of in- 
cluding the H 3 terms in the equations of motion. 

The first step is to find the transformation that uncouples H 2 . This can be done by regarding 
as a matrix, 


l 

1 q 


:s] 
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where 


and 


q ~ < 


<*3 

Pi 

p 2 


y , a column vector, 


IS] 


l " T 7 
F 2 
0 

0 

0 

1 “ r 0 


0 

1 +k -kr 
0 

1 

0 

0 


0 

0 

k ( r 2 ~ r i) 

0 

0 

0 


1 - r. 


1 0 
0 0 
1 


0 

0 


1 



0 


0 

0 

0 

0 

0 

1 


and finding a constant -coefficient matrix [D] such that, for q", a column vector in the new 
coordinates, 


q - [D] q" 

represents a canonical transformation to a system H** (q", p", t), for which H** is a diagonal 
matrix. This is done in two steps, first a transformation matrix [K] to a system H* (q' , p' , t ) in 
which 


3 



i = 1 


and then a transformation [m] to a system H** (q", p", t) in which 


H 


2 


\ r. ( p i') a + ( q i ") 2 • 


Then the desired matrix [D] is 


Ed] = [K] m . 
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The procedures used are quite general and could be applied to any similar problem. In the par- 
ticular case under consideration, the terms p 3 and q 3 in h 2 (q, p, t) are already in the desired 
form. Thus a matrix [D'j which diagonalizes the matrix [s'] will be developed where 



The prime notation will be dropped in the remainder of this section; however the meaning is 
obvious. 

A formal procedure for obtaining the first canonical transformation is developed by Pars 
(Reference 29); only the essential steps are outlined herein. This part of the problem can be re- 
stated as: Given a real symmetric matrix [S], find a matrix [K] such that 


[K] T [S] [K] 



where [L] is a diagonal matrix 



and such that 


[K] T [Z] [Kl = [Z] 


where 

0 10 

0 0 1 

0 0 0 

-10 0 

Then the matrix Ck] will define a contact transform 

q “ [K] q 
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which accomplishes the desired result. The first step in obtaining the matrix [k] is to find the 
eigenvalues of the matrix [z] Csl (i.e. find the roots k of Is +^z| - 0). 

This is equivalent to finding the values of k for which the determinant equation 


1 - r. 


0 

k 

1 -r„ 


0 

1 +k -kr 

1 




1 - r, 

T ~2~ 

-k 


r~ o 


= o 


is satisfied. This equation can be expanded to give 


A . 4 + \ 2 


p-r 2 14 

L'i' a ' + 


k - kr, 


- 2 


(i' r 2 ) 


1 + k , . . 

+ r "r 2 (l- r i) (l' r 2 ) 


(4-4) 


(4-5) 


The roots of this biquadratic equation are the normal frequencies of the two decoupled normal modes 
of the linear problem associated with h 2 . When e = 0, k =3 and Equation 4-5 can be reduced to 


r 1 r^ ( U2[-3r 1 U 2 r 1 r ! *2 ri -r jt l]H(l-r 1 )(l- r! ) = 0. 

which (with the substitution k - i-v) is the same as Equation 4-4 of Breakwell and Pringle. The 
roots of this equation occur in complex conjugate pairs; when 0 < r x < r 2 <l, both the roots of 4-5 
and the roots ±\ 3 associated with the characteristic equation for the already decoupled pitch motion 


A 2 + k^r 2 -r 1 ) - 0 (4-6) 

are pure imaginary numbers. When k has any of the four values +iw 1? +ia> 2? -i co l , -i^ 2 (where eo 2 ><, 1 ) 
obtained from Equation 4-5, there is a nontrivial solution to the matrix equation 

[^{cj = i = 1. ■•',4 

where the C i are the eigenvectors of [zl [s] . A set of eigenvectors for the matrix can be obtained 
from the minors of any column of the associated determinant; see Wiley (Reference 30). One set of 
eigenvectors obtained from the second column of Equation 4-4 is based on 

M 1 _r i ~ r 2 ) 

r, - 1 + r 2 
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(l + X 2 ) {r 1 -r 1 r 2 ) 

X [-^ 2 r i r 2 ‘ ( 1_r 2 ) ( 1 _ r l )] 


and in particular is given by 


c i = 




iw i l 1 = r i _r 2 ) 

r i " 1 + r 2 

(1 -^i 2 ) (rj-r, r 2 ) 
[a- 2r i r 2 - (l-r 2 ) (1 - r i )] 

-r 2 ) 

cm i 2 r 1 - 1 + r 2 

(i - &*) ( r i - r i r 2 ) 

h ia i [ &, i 2r i r 2 _ ( 1_r 2) ( 1_r i)] 


c, = cr< 




iw 2 (l" r i ~ r 2 ) 

“ 2 2 r i _ 1 + r 2 

(i -« 2 2 ) (r t - r j r 2 ) 
ia; 2 r t r 2 - (1 - r 2 ) (l - r^J 

-ia 2 (1 - r x - r 2 ) 
r l - 1 + r 2 

( 1 _<iJ 2 2 ) ( r i _r i r 2 ) 

-io .' 2 jo .' 2 r t r 2 - (l - r 2 ) (l - rjj 


where p, o, y, and 8 can be any arbitrary (real or complex) constants. Now define K as the matrix 
formed from the columns C 15 C 2 , C 3 , C 4 , in that order. Then 


M T Lz\ [K] = [F] , 

where 

0 

f 2 
0 

0 


[Fl 


0 

0 

0 

-f„ 


f, ■ {£,} t la {£,} . 


f , 1 Kt 1 [zl { c 4 ■ 
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Now, if the values of p, a, y, and s are chosen such that f x = f 2 = 1 , then, in 

M T [Z] [K] = [Zl , 


the matrix [K] is symplectic and the transformation 

q = [K] q' 


is a contact transformation. In addition, 


[K] t [s] [K] = [E] , 


where 


and thus 


[E] 


H 2 = J {q} T CS] {q} 

= y{q'} T [K] T [S] [K] {5'} 
= ■§• {q'} [E] {q'} 


- \ qj Pi + A. 2 q 2 ' P 2 

Thus the matrix [k] defines a contact transformation which transforms H 2 to H* and satisfies 
Equation 4-6. The requirement that f 1 = f 2 = i provides two relationships for determining p, o-, 
y , and 5 : 

2Si "2 {(iTj -r 2 ) (l-^, 2 ) (r t - ri r 2 ) + (l- r 2 - r 2 - (1- r 2 )(l - ri )]| = 1 , (4-6) 
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and 


2^Sio; 2 |(l-r 1 -r 2 ) (l-^ 2 2 ) (rj-fj r 2 ) + (l-r^Fj^a) r 2 - (l-r 2 ) (l- ri )]| - 1. (4-7) 

Any values of p, a, y , and 5 that satisfy the above equations are satisfactory in defining [K] . Thus 
it would be possible to make two arbitrary choices, such as p - 1, cr = 1, and find values of 7 and 
5. However, it is more convenient to obtain the second transformation matrix [M] and then find 
two additional relationships that simplify the final transformation matrix [D] . 

The next step in the process is to find a canonical transformation of the form 

{?'} = [M] {q"} 

that transforms h 2 *, (q.', p.' ) into h 2 **, (q." , p.")? where 


and 


(q'. p') 


2 



i = 1 


H**(q",p") " £pi" !+ “i 2 q; ! ' 

i = 1 

There is a generating function W, given by Whittaker (Reference 31), which can be used to find the 
needed transformation. Consider 



i = 1 


with 




gw 


1, 2 , 


and 


Pi' 


<9W 


dq/ 


i “ 1,2. 
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I 

i 

I 


The above leads directly to the relations 

„ Pi" 

+ X~ * 

i 

, *•! q i" £ L 

p i “ 2 + 2 ’ 

which can be expressed in matrix form as 


{q'} - [M] {q"} , 


where 


1 0 

o 1 

\ 

2 0 

\ 

0 “ T 


The first transformation matrix [K] is of the form 


[K] 


PC 1 1 

oC 21 

~yCi x 

- SC 2 ! 

PC 1 2 

°C 22 

yc t 2 

sc 22 

P^13 

° C 23 

^ 1 3 

SC 23 

PC 14 

aC 24 

- 7 C 14 

~ SC 24 


and thus the matrix [d] would have the form 


[D] = 


p + ~2/ C n 

Vy 

p 2 / w i 2 


p 2/^13 

P + 2 


k 2 B 

+ "2 ) C 2l 


i a 


a + 


M c 

2 / C 23 
2 ) C 24 


(< 4 

M- 

(M 


14 


M' 

(£4 

M) 
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The freedom which exists in choosing p, <?, y, and S permits the choice 


_ 

y - K 




which simplifies [D] to the form 


[D] 


2pC x l 2 ctC 2 t 


0 

0 

2 pC 


0 
0 

14 2ffC 24 


0 

2/QC 12 

0 


0 

2crC 


22 


s 

2oC 


23 


0 


and Equations 4-6 and 4-7 become 

4p 2 {(l -r, -r 2 ) (l -to*) (r, -r x r 2 ) + (l - r 2 ~ r^ 2 ) [u 2 r t r 2 - (l- r a ) (l ~ r , )] J- = 1 


4p 


! {( X - F i ' r 2 ) ( 1_ "2 2 ) ( r i " r i f 2 ) + [^2 2 r Z r 2 ~ t 1 - r 2 ) ( X _ r i)]} = 1 


(4-8) 

(4-9) 


These expressions are rather unwieldy; however, they can be reduced with some algebraic sleight- 
of-hand. Write Equation 4-8 as 

4p 2 ^(l-r 2 -r 1 o. 1 2 ) + {(l " r 2 - r, a 2 ) r^ 2 - (l - r 2 ) (l-r,)] 

+ (l -t 1 ~r 2 ) (l -aij 2 ) (rj-r, r 2 )}'] = 1 . 


The term inside the brace can be written as 


-i -, 2 r 2 ^ 2 a , 2 2 + 


( 1-r 2 ) r i r 2 K !+ “ 2 J ) + ( 1_r 2 ) [- r i 2 - 2r i r 2 + 2r 1 - 1 + r 2 ] , 


(4-10) 


and a> 2 and a>* are the roots of this equation, which can be written as 


r i r 2“ 4 + 


, „ 0 foe 2 45e 4 \ , 

3r l 2r i r 2 ~2 r i - i + r 2 - [ 2 + § / r, (l -fj 


4 ( 1_r i) ( 1 “ r 2 ) + (lT + ^“) ( 1_r 1 ) (l- r 2 ) " °- 
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Thus 


CO 


2 

1 


+ 


CO 


2 = 
2 



r 2- 2r i 


1 + r _ 



r 


l 


r 


2 



and 


4(1- 



r i) ( 1_r 2 ) 


When these expressions are substituted in expression 4-10, it becomes identically zero; thus 
Equation 4-8 can be written as 


4 P 2 (l - r 2 


r i a i 2 ) r i r 2^i 2 ~ r ‘i) 


i , 


and by a similar process Equation 4-9 becomes 


4cr 2 (l - r 2 - r t c a 2 ) r t r 2 (co* - ^, 2 ) 


1 . 


These two equations can be solved to give 


2 P 



i 



and 


2 cr 


[(! - r 2 


i 



With these two substitutions, [D] becomes 
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&, l( r i +r 2" 1 ) 

^( r i +r 2 _1 ) 



[( 1 " r 2" r i"l 2 ) r i r 2K 2 " a, 1 2 )] 1/2 | 

_( 1_r 2~ r i"2 2 ) r i r 2K 2-a, 2 2 )] 1/2 

0 

0 

0 

0 

-(l-r 2 -i-, 2 ) 1/2 

-(l-r2-r x ^)^ 

“l^l^K""! 2 )] 172 

° J 2 Pi F 2 ( M 2 ~ 01 1 )] ' 

0 

o 

(i-^ 2 ) r i (l - r 2 ) 

(l-^)r 1 (l-r 2 ) 

«,[l 

[ 1_r 2- r l a i 2 ) r ir 2 K""l 2 )] 1/2 « 2 [( 

1_r 2- r i"2 2 ) f i r 2h 2 --2 2 )] 1/2 

a 'l[( 1 " r 2)( 1_r i)" r i r 2"l 2 ] 

W 2[( 1 “ r 2)( 1 “ r i)- r i r 2 W 2 2 ] 



[( 1_r 2 _r i a -'l 2 ) r i r 2K 2_a)2 )] 1/2 

[(! - r 2 - r i ~ 2 2 ) r ! r 2 K 2 -" 2 2 )] 1/2 

0 

0 


( 4 - 11 ) 



Thus the matrix Dd] simplifies the roll-yaw portion of the Hamiltonian. 

It is also possible to find new variables that replace q 2 and p 3 , and eliminate the H 0 and H 1 
terms. This time, consider the following portion of the original Hamiltonian: 

H' = T t (t) - 1 + p 3 Tj (t) + (p 3 2 + <x> 3 2 q 3 2 ) . 

Two successive canonical transformations, the first derived from a generating function: 


F ( C >3’ P 3 ' ’ *) " ' b ( 1 > P 3* ■ 

with 


dF_ 

p 3 dq p 3 ’ 

, _ d¥ 

q - r - q - b( t ) , 

*p 3 


and 


H* - H' -bp # , 


and the second derived from a generating function 


with 


F *(q 3 ". p 3 '. t) - -p 3 ' q 3 " + C(t)q" 


j" jbc — (c 2 + <u , 2 b 2 ) - Tj ( t ) + lj dt 


dF 
*P 3 ' 


- - q 


= - il. 

dq" 


- P,' - c(t) 
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and 


= H* + bq" + 


j^bc - y (c 2 + a > 3 2 b 2 ) “T x (t)+lJ , 


combine to give a transformation 


q 3 - q 3 " + b(t) , 
p 3 3 p 3 " + c ( t ) • 


and a new Hamiltonian 

H ** = T (p 3 " 2 + “ 2 2 ) + p 3 " [° ' + T i (t )] + q 3 " K 2 b + i ) ■ 

Now if c(t) and b(t) are chosen so that 


c ( t ) - b(t) + Tj (t) = 0, 


and 


c( t ) + co 2 b( t ) - 0 , 


i.e., if c(t) and b(t) are solutions to 


c(t)+a- 2 c(t) = - ^ 3 2 T i ( t ) , 


and 


b(t) + a; 2 b(t) - T 1 (t) , 


then the Hamiltonian H** becomes 


(q 3 " . p 3 " . t) - y ( p 3 : 


+ <^>-1 


with Tj (t) defined by Equation 4-1, b(t)becomes 


b(t) 



s in ( it ) , 


(4-12) 
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where 


2e e 3 

1 " co 2 -l 4(co 2 -l) ’ 

5e 2 _ lie 4 

2 " co 2 - 4 6 (co 2 -4) ’ 

39e 3 

3 ~ 4 (co 2 - g) ’ 

103e 4 

4 6(co 2 -16) 

Similarly 


c(t) 

where 

2e e 3 co 3 2 

Cl ' co 3 3 -1 " 4(co 3 2 -l) ' 

5e 2 lle 4 <^ 3 2 

C2 " 2(co/-4) " 12 (co 3 2 -4) ’ 

1 3e 3 co 3 2 

C3 " 4(co 3 2 -9) ’ 

103e 4 co 2 

" 4 = 24(o> 3 2 - 16) 

When the above transformations are all applied to the original Hamiltonian given by Equation 4-3, 
the new Hamiltonian has the form 


L 


c . cos ( it ) 


(4-13) 


H**(q\p",t) = 2_V*(q\p",t: 
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where 


H 0 **(q", p", t) = 0, 

H**(q",p",t) = 0, 

3 

H**(q",p",t) = (Pj" 2 + w j 2 qj" 2 ) ■ 

i - 1 

and H 3 ** is obtained by making the following substitutions in h 3 (q, p, t): 

= d u q x " + d 12 q 2 " > 

^2 = d 23Pi" + d 24P 2 " ■ 

P 1 ” d 3 3 P 1 T d 34 P 2 * 

P 2 ” d 4 1 ^1 + d 4 2 ^2 * 

< 3 3 = ^3 + b (*) • 

P 3 * P 3 " +C(t), 


where d ij represents the element in the i th row and the j th column of the matrix [D] defined by 
Equation 4-11, and b(t), c(t)are defined by Equations 4-12 and 4-13. Thus these new coordinates 
lead to a system which has no terms of order 0 or 1. The terms of order 2 will be removed, in 
the next section, by a final canonical transformation. 

Solution of Second Order Equations 

In the preceding section the original Hamiltonian was considerably simplified by a series of 
canonical transformations. The end result was a system in which the linearized equations repre- 
sent three uncoupled harmonic oscillators. The solution to the linear equation is well known; how- 
ever, when it is obtained by finding Hamilton’s principal function S, then this function generates a 
contact transformation which will remove the second -order terms from the new Hamiltonian. 

The process of finding the principal function for a Hamiltonian of the form 

H ** = £ ![(».”) 
i = 1 
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follows directly from the solution for a single harmonic oscillator (Goldstein, Reference 32). The 
principal function S(q", a, t) is 


= L{{h-^' 2 ] dq ‘ " a ‘ *} 


and the transformation defined by this function is obtained from 


a = 

P i da. 


a) i 


— arc sin \ t 


which is equivalent to 


and 


J2a[ 

di" = sino>. (t +/3.) 


Pi 


dS 

d *i " 


y2a~ { cos a), (t +/3.) 


The new Hamiltonian, H***, is related to the old one by 


H** + _|S_ = H ** _ £ 0 . 


Thus the new Hamiltonian is 


(a, {3, t) 


■ IX 


(a, 0, t) 


(4-14) 


where 


Hq*** ( a , Pi t ) = 0, 

H p, t) = 0, 

H***(a, P, t) = 0, 
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and H*** (a, /?, t) is obtained from Equation 4-3 by making the following substitutions in h 3 (q, p, t) 


V2 d tl 


Z 5 7 sin "i (t + ^i) + 


i2 d 12 


(* + ^ 2 ) 


q 2 - ft d 23 -fTl cos O.J (t +/3j) + /2 d 24 i °-~2 cos ^2 l* 1 + ' fi 2 ) 
Pi = ft d 33 ftl COS “l (t + ^l) + ft d 34 /oj COS C 2 (t +/3 2 ) 


^ d 41 


V2 d 4 


j i^T sina i (t + ^i) + u - 2 ftl sin&: 2 ( 4 + ^ 2 ) 


V2 


"*3 6 c' 


“V 5 ^ sina 3 (t + ' fl a) + C ( t ) - 


ft ft~3 COS a 3 (t + ' 3 s) + bC 4 ) ■ 


(4-15) 


Summary 

In this chapter, the equations describing the attitude motion of a rigid gravity -gradient- 
stabilized spacecraft moving in a known eccentric orbit has been formulated from the Hamiltonian 
point of view. A series of canonical transformations, following Breakwell and Pringle (Reference 
12) were developed; these considerably simplify the Hamiltonian. The formulation is extended to 
include terms to fourth-order in eccentricity. 

The final form of the Hamiltonian given by Equation 4-14 can be easily interpreted. When 
only the linear terms are included (i. e., when H*** and all higher -order terms are deleted) the 
Hamiltonian equations 


dH*** (a, t) 

da . 


(4-16) 


and 


^H*** (a, j°, t) 


have the solutions 


P, (t) = P { (0) 
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and 


a. (t) = a. (0) . 

The original variables of the problem (i.e., the and p £ ) are then given by Equations 4-15 
as combinations of sinusoidal oscillations with constant amplitude and frequency. In fact, the 
amplitudes are determined by the initial conditions, and the frequencies are the eigen-frequencies 
or frequencies of the normal modes found from Equations 4-5 and 4-6. Unbounded motion only 
occurs when oj 3 = 1, for in this case b(t) and c(t) are unbounded, as can be seen from Equations 
4-12 and 4-13. 

The situation changes when the nonlinear terms are included. The third-order part of the 
Hamiltonian, h 3 *** (a, /3, t) is not identically zero; therefore Equations 4-16 no longer have a 
constant value as a solution. The terms in these equations are small in magnitude, either 
because they are second-order in small displacements and momenta; because they involve a 
small parameter e ; or for both reasons. Thus the motion may be regarded as consisting of 
sinusoidal oscillations with slowly varying amplitude and frequency. In the next chapter the 
method of averaging is applied to find regions where these small variations will have a signifi- 
cant effect. 
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CHAPTER 5 


i 


GENERAL PARAMETRIC AND NONLINEAR RESONANCES 


Introduction 

An approximate solution to the single-degree-of-freedom, linear, in-plane, pitch equation was 
developed in terms of asymptotic expansions in powers of a small parameter in Chapter 3. A sec- 
ond form of this general approach, called the method of averaging (or of Krylov -Bogoliubov) is 
applicable to multi-degree-of-freedom equations such as those of the previous chapter. The theory 
for both these methods of solution is given by Bogoliubov and Mitropolsky (Reference 27). A par- 
ticular implementation of the method of averaging is used herein; it follows closely that of Breakwell 
and Pringle (Reference 12). The formulation has been extended to include all terms of second or 
third order in the generalized coordinates and through fourth order in eccentricity. 

The approach outlined in Reference 27 is to consider the problem as defined by kinetic and 
potential energies (which are quadratic forms in the generalized coordinates) and to assign the 
nonlinearities to a perturbation term. This implies that both the nonlinear effects and the other 
components of the perturbing term are "small.” The equations are transformed, first to normal 
coordinates and then again to reflect the harmonic solution of the linear portion of the problem. 

Once this is accomplished, the equations are in what Bogoliubov and Mitropolsky refer to as "the 
standard form," which is 

d x k 

~3t = £ \{ t - x r x 2 - x n) - k = 1. ••■.n, (5-1) 

where e is a small parameter and X k may be represented by the sum 

X k (t. x !. •••-* X„) = eXP ( ivt)X kv( X l- " ' ’ X n) ' k = 


This is exactly the result of Chapter 4 in which the Hamiltonian is developed through third 
order in the generalized coordinates and momenta, and through fourth order in eccentricity. The 
transformations to normal coordinates q", p", and then to coordinates a, /3 (which reflect the simple 
harmonic motion associated with the unperturbed linear equations) are accomplished through the 
applications of canonical transformations. 

In the specific case under discussion. Equation 5-1 is replaced by the six Equations 4-16 de- 
veloped in the preceding chapter. The intent of this chapter is to find, by the method of averaging, 
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the portions of the Lagrange region where the perturbing terms due to orbital eccentricity and 
nonlinear coupling have a significant effect on the motion. 

The basic assumption of the principle of averaging is that the coordinates of the problem un- 
dergo a motion composed of a slowly varying term and small vibrational terms. That is, the time 
required for any of the variables to change appreciably is long compared with the oscillation period 
of the unperturbed coordinate. When this is the case, Equation 5-1 is replaced in the first approxi- 
mation by an equation of the form. 

r , x-i 

~dt ~ eI ^L X k( x i’ x 2’ ***’ x Jj ■ k = 1, 


where 


x k - + small vibrational terms , 


M is the operator of averaging with respect to time, and the f are held constant during the averag- 

t 

ing. The theory can be continued to develop higher-order oscillations that reflect the effect of the 
small oscillations in x on the perturbing force and thus on the motion in general. 

There are no constant terms in h 3 *** (a, j3, t) ; thus, when the Hamiltonian is averaged over time, 
its value is zero and the coordinates represent unperturbed motion except in a limited number of 
cases. The third- order Hamiltonian contains a number of terms that are themselves products of 
harmonic terms. For example, the term 1 - r 2 /r 2 qi p 3 p 2 that appears in h 3 (q, p, t) is transformed 
into 


1 7^ [ d H d 41 K') 2 + ( d ll d 4 2 +d 12 d 4l) +d !2 d 42 KT] [ P 3 ' + «<*)] 


in h 3 ** (q\ p", t) and then into 


~ T V' [ 'l/' 2a i sin2 "x ( 4 + / 3 l) + ( d n d 42 + d 12 d 4l)73^ V a l a 2 sin<i, i (t+^Jsinco^t +/3 2 ) 

( 4 + / e 2)][/2°7 cos a) 3 (t + £ 3 )+c(t)J 


d 12 d 42„ -2 

+ 2 — 2 a 2 sir\ z 0 ) 2 


in h 3 *** (a, /3, t ) . There are two types of terms in an expression such as that given above which 
may not integrate to zero. They become more apparent if it is noted that 

sin 2 ^(t +/3 1 )cos<u 3 (t + p 3 ) = cos o) 3 ( t + /3 3 ) 

- ^ jcos [(2wj +*J 3 )t + +co 3 /? 3 ] + cos [(2^ 1 -co 3 )t+2<u 1 ^ 1 -<2) 2 /3 3 ]| - 
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and that 


sin 2 (t +/3j) c(t) 

“ ^ { cos [( 2aJ i + i)t + 2a. j /jJ + cos [( 20 .J - i) + 2^ /?j| • 

i= 1 

Thus the first term does not average to zero when 2co 1 = oj 3 and the second term does not 
average to zero when 2 ^ = 1, 2, 3, or 4. These are examples of internal and external resonances 
respectively. This nomenclature arises because the internal resonances are inherent in the sys- 
tem, even for a circular orbit. They originate in the expressions for the kinetic and potential ener- 
gies, while the external resonances arise from commensurability of frequencies of the external 
forces and accelerations with the normal frequencies of the system. When the orbit is circular 
(i.e., when e = 0), all the coefficients a. , b., and c. are zero and there are no external resonances. 

Breakwell and Pringle demonstrated the existence of two internal and four external resonances 
for spacecraft with inertia parameters in the Lagrange region. The internal resonances appeared 
for a> 3 % 2 o) 1 andaj 3 % OJ 2 ~ co i • External resonances were found when a> 1 or ^ 3 = 1/2, when o> 3 = 1, 
and when a> 2 -co x - 1/2. A number of additional external resonances become apparent when terms 
through the fourth order in e are included in the expansions for a(t), b(t), and c 2 (t). Additional 
resonances, both external and internal, could also be obtained by inclusion of fourth-order terms 
in the Hamiltonian, or by consideration of higher-order approximations in the method of averaging. 
The algebra involved in either of the latter two steps would be almost prohibitive, and they are not 
included herein. 

The method of averaging presents a set of first-order differential equations whose solutions 
define, within some degree of approximation, the time variation of the parameters (a., /3.) and 
thus the original variables (q. , p.) . However, it is not necessary to obtain the solution to these 
equations; also, the approximations in the analysis might cause misleading results. The analysis, 
instead, is intended to determine to what extent the variations are bounded and what conditions are 
required to excite significant changes in the unperturbed motion. 

Resonant Frequencies 

Any rigid gravity- gradient spacecraft with inertia parameters in the Lagrange region can be 
shown to be Liapunov- stable for a circular orbit, as shown by Pringle (Reference 16). However, 
there are inertia parameters for which the averaged third-order Hamiltonian is not zero. These 
statements do not conflict; for the Lagrange region, internal resonances have been shown (Ref- 
erence 12) to produce bounded, but occasionally large, interchanges of energy between pitch, roll, 
and yaw motions. The situation changes markedly when the external resonances are considered 
for spacecraft in the Lagrange region. In these cases orbital eccentricity, however small, can 


4 


c . cos ( it ) 
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excite unbounded oscillations in one or more of the spacecraft modes. This is just a three- 
dimensional generalization of the results obtained in Chapters 2 and 3 for the one-dimensional 
pitch-motion case. 

In discussing the situations in which resonances occur, it is helpful to remember that the a>. 
are the frequencies of the normal modes of the linearized equations of motion. They are obtained 
from Equations 4-5 and 4-6, which can be written with k = ico as 

f l r 2 a ’*2 + [ kr i 2_ 2r i r 2 "( k ' X > r i + r 2 _ + ( k+ ^ (* _r i) ( 1_ T l) = 0 


and 


“i ■ k ( r 2 _r i) = 0 • 
Thus, since by definition w 1 <co 2 , therefore 


[- b - /b 2 - 43c"! 1 

= [■ 21 J 

T- b + /b 2 ' 4ac 

L J 


(5-2) 


(5-3) 


and 

“a = [ k ( r 2 “ r i )] I/2 > ( 5 “ 4 ) 

where 


b - k r ^ “ 2rj r 2 - (k~ 1) r x + r 2 - 1 , 
c = (k + 1) (l- r t ) (l - r 2 ) , 


and 


k 


3 1 + 


3e 2 

2 


+ 


15e 4 

8 


)• 


When the Hamiltonian h 3 *** (a, /?, t) is expanded completely in the manner illustrated above, 
the combinations of frequencies shown in Table 5-1 are resonant. 
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Table 5-1 


Resonant Frequency Combinations. 




Resonance 



Effect 


Resonance 



Effect 

1. 

"3 " 

1 = 0 



D 

17. 

" 3 + 2 "l 

- 2 = 

0 


D 

2. 

2 co l 

-1 = 0 



D 

18. 

" 3 " 2c °2 

+ 2 = 

0 


I 

3. 

2oj 3 ■ 

-1=0 



D 

19. 

"3 + "l 

+ oo 2 - 2 

= 

0 

D 

4. 

G) 2 

"i - 1 = 

0 


N 

20. 

W 3 + “2 

- a> l - 2 

= 

0 

N 

5. 

“3 - 

2 cjl> x = 0 



I 

21. 

oj 3 - o» 2 

- oo 1 + 2 

= 

0 

I 

6. 

"3 + 

“l “ “2 = 

0 


N 

22. 

2 oo 2 - 3 

= 0 



D 

7. 

"3 " 

2 o> 1 + 1 = 

0 


I 

23. 

2" 3 “ 3 

= 0 



D 

8. 

"3 " 

2co 1 - X - 

0 


I 

24. 

"l + “2 

- 3 = 



D 

9. 

"3 + 

2co 1 - 1 = 

0 


D 

25. 

oj 3 + 2o>, 

- 3 = 

0 


D 

10. 

"3 + 

CO 1 - OJ 2 + 1 

= 

0 

N 

26. 

*>. 3 + 2C0 2 

- 3 = 

0 


D 

11. 

“3 + 

"2 - "l - 1 

= 

0 

N 

I 

27. 

" Za} 2 

+ 3 = 

0 


I 

12. 

W 3 

CO 2 - co 1 + 1 

= 

0 

i 

I 

! 

28. 

W 3 + “l 

+ co 2 - 3 

= 

0 

D 

13. 

2<Wj - 

-2 = 0 



D 

29. 

^3 + ^2 

— co^ — 3 

= 

0 

N 

14. 

2o> 2 - 

-2 = 0 



D 

30. 

2o) 2 - 4 

= 0 



D 

15. 

2w 3 - 

-2 = 0 



D 

31. 

2 °> 2 + “3 

- 4 = 

0 


' D 

16. 

"l + 

CO 2 - 2 = ' 

0 


D 

32. 

“l + “2 

+ o> 3 - 4 

a 

0 

i 

D 


D indicates divergent oscillations result 

I indicates an interchange of energy between the pitch and roll -yaw modes 
N indicates no significant effect 


The first six terms have been discussed in Reference 12 and with the exception of the linear 
pitch resonance case (co 3 = 1) are all either second-order in the coordinates and first-order in ec- 
centricity or third-order in the coordinates. The next ten items are third-order in the coordinates 
and first-order in eccentricity or second-order in both coordinates and momenta. Items 17 through 
24 are either third-order in eccentricity and second-order in the coordinates or the reverse. The 
resonances for items 25 through 30 are third-order in both coordinates and eccentricity or second- 
order in the coordinates and fourth-order in eccentricity. The resonances for items 31 and 32 are 
both third-order in the coordinates and fourth-order in eccentricity. 
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Each of these resonant-frequency combinations defines a line or point in the inertia space of 
the Lagrange region. These lines are plotted from numerical evaluation of Equations 5-2, 5-3, and 
5-4 with k = 3 (i.e., for e = 0). The first six lines are shown in Figure 5-1; the next ten are shown 
in Figure 5-2, the following eight in Figure 5-3 and the last eight in Figure 5-4. The same combin- 
ation of lines is shown in Figures 5-5 through 5-8 respectively where the equations are evaluated 
with k = 3.189 (i.e., for e = 0.2). The similarity of the figures demonstrates the relative insen- 
sitivity of the locations of most of the lines of resonance to variations in e. 

Behavior Near Resonances 

The Hamiltonian H*** (a, /3, t ) has terms with a steady non-zero component whenever any of 
the 32 resonant relationships of the preceding section are satisfied. Whenever this occurs, the 
averaged perturbation equations 


°i = - 2^7 (a. /S, t )] 

and 

: 2^7 M[H*** (a, /3, t)] 


are not identically zero; thus, there are long-term changes in the values of the a. and /3. . It is pos- 
sible to determine the nature of these long-term variations without obtaining formal solutions for 

a. (t) and /3. (t). 

The procedure is best explained with the aid of several examples. Consider the case where 
2a x - l r £ ; then the slowly varying terms in the Hamiltonian are 


(a, /3, t)] 


d 1 1 ( d 1 1 + d 4 1 ) 


+ d 23 ( d 2 3 + d 3 3 ) I 


C 1 + 


3(1- rj 


d 2 s 
U 23^1 


+ 


3(1 


r 2 ) d ll ^23 


[ b 1 + \ ( S 1 b 2 



l)t + 2a.. /3j] 


(5-7) 


Each of the terms b 1? c,, and s t , are first-order or higher in e; thus Equation 5-7 can be written. 


H***(a,/3, t) - eKa-j cos (et + 2a t /3 t ) . 


It is possible to obtain the same equations for Oj and from a Hamiltonian that is independent of 
time and thus a constant of the motion, by considering new variables a*, p* such that 


»i* (t) = «i(t) 
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Figure 5-3— Resonant-frequency combinations 17-24 in Table 5-1 (e = 0). 
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Figure 5-8— Resonant-frequency combinations 25-32 in Table 5-1 (e = 0.2). 
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and 


/vm = Pi <t> + ■£[ 

together with a new Hamiltonian h' (a*, /3‘) defined as 


H'(a*,/3*) = eKa; cos + 2 ^ • 


The perturbation equations are 


a* - 2^ 1 eKa* sin 2^^ 


and 


Z 3 !* 


eK cos 2^j /3 t 


e 

2^; 


The nature of the solution to these equations depends on the magnitude of e relative to e. When 
I e/ 2 cv 1 1 > | eK | , j3* is always positive for e > 0 and negative for 6 < 0; thus p* (t) increases monot- 
onically for e > 0 and decreases monotonically for e < 0. In this case a* (t) alternates in sign and 
a* (t) exhibits bounded oscillation. 


When | eK | > 



lim /3 * (t) 

t — 00 x 


arc cos 


e 

2eu^K 


in such a manner that lim a * (t) -*a>. This can be more easily seen from the phase plane considera- 

t -CO 1 

tions which will be discussed next. 

The new Hamiltonian h' (a * , p * ) is a constant of the motion, thus it is possible to draw lines of 
constant H in a phase space with coordinates a* and 2^ x p* . This is done in Figure 5-9 for the case 
| eK | < \e/2cv\ and in Figure 5-10 for | eK | > |c/2wj. Figure 5-9 is drawn for e > 0 and Figure 5-10 for 
K > 0; it is only necessary to change the direction of arrowheads for the cases e < 0 and K < 0. 

There is a third approach to the problem which also depends on the constancy of H' (a*, p *) . 
Equating values of h at t = 0 and t = t gives 


cos 2co 1 P* (t) 


H 0 ' ^ « i* ( t ) 

eK 2e co t K 

(tj 
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where h 0 ' = H' (a*, £*) evaluated at t = 0. However, a x (t) is zero only when cos 2 '& X P X = ±1 (i.e., 
when sin 2 o3 1 /3 1 = 0). Thus (t) = 0 only when 

# H o €a i 

±a i eK 26^^ 


The two sides of this equation are plotted versus a* (as firm and dotted lines, respectively) in 
Figure 5-11 for e/ 2 e a> 1 K > 1, which clearly reveals the extent to which a* (t ) is bounded. Fig- 
ure 5-12 makes the same plots for 0 < e/2e a> x k < 1. In this case it is apparent that if a* (t) is ever 
positive it is always positive, and it can be seen from Equations 5-9 and 5-10 that a* (t) will be- 
come positive for any initial conditions. 

Each of the three approaches leads to the same set of conclusions. The value of a * (t) in- 
creases without bound for 2*> 1 sufficiently close to 1, and the minimum frequency mismatch (i.e. e) 
required to guarantee bounded motion increases linearly with eccentricity. As has been noted be- 
fore, the first conclusion means that the motion increases until higher-order terms become 
important. 

Fortunately, this same analysis applies to a number of the other resonances shown in Table 5-1. 
In particular, when 2a 3 “lie, the slowly varying terms in the Hamiltonian have the form 

(a, /3, t)] = eK' a 3 cos [^( 2 oj 3 - l) t + 2^3 /3 3 J , 


where k' can be obtained as a collection of terms in the same manner as was done above. Obvi- 
ously the same conclusions apply to this case as did for the previous one. In fact, when any of the 
frequencies are such that 2cv. - j ¥ e , the Hamiltonian will have the form 

m[h*** (a, /3, t)] = e J Ka. cos [(2a:. - j)t + 2co. /3.J , 

where K may again be obtained by collecting the appropriate terms in the expansion of H***(a, /3, t), 
and the width of the resonant band is proportional to e j . Thus frequency combinations 13, 14, 15, 22, 
23, and 30 all lead to unbounded motion if the appropriate value of e is sufficiently small. This is 
the three-dimensional generalization of the parametric resonance phenomena studied in Chapters 2 
and 3. These unbounded notions arise from the variation in the restoring torque with orbit position, 
rather than from the nonlinear terms. 

Both the combinations cc l +oj 2 = 2, and +co> 2 = 3, have long-term variations of the form 

M[H(ol, J3 t t)] = Ke j cos [(^ + “ j j 1 +0J i @i + 03 2 ^ 2 ] ' 
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-Although this form differs from any that have been previously discussed, it is similarly applied. 
With ojj + co 2 ~ j — e, the averaged Hamiltonian can be written as 

H*** (a.,/3, t) = Ke j a 2 cos [<^ 2 /3 2 + oj 1 J3 1 + € t] . 

The pitch variable <x 3 is not included in the Hamiltonian, indicating that, except for small oscilla- 
tions, it is constant. A new Hamiltonian which is also a constant of the motion can be developed 
by considering new variables 


a 


l 


a 


1 ’ 


and 


+ a > 1 


and a new Hamiltonian, 


H' = Kei cos [o> 2 P 2 + *>j £*] + — 

The canonical equations for a* (t) and a 2 ( t ) are 


^Ke 1 


V“ 1* a 2 Sin [" 2^2 + "l ^1 ] 


and 


a 2 - co 2 Ke } ^a*a 2 sin [^ 2 + a> l J3* ] ; 


thus, 



and 


a 


2 


( a i* + c i) 
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The new Hamiltonian can be written with the above as 


€a x 

V = Kei y — fa; (a* +C) cos [a> 2 £ 2 p*] + ■ 


and a* (t) is only zero when 


H o- ^7 

c = — 

Ke> 

A graphical construction of the two sides of this equation would be quite similar to Figures 5-11 
and 5-12, except that f (a*) = j ** 2 + a *c x instead of f(a 1 *) = a * . Thus, for any initial conditions 
such that a* (0), a 2 (0) ^ o, a* (t) is unbounded if | e/V ~ 2 Ke j | < 1. Since 





+ a. 


2 2 1 
“2 (t) = ^ tt l (O + ’ 

this means that both the roll-yaw modes are unbounded, to the order of approximation of this so- 
lution, when e is sufficiently small. 

The 19th, 28th, and 32nd combinations of frequencies in Table 5-1 are all of the form 

<^1 + ^2 + ~ j = 0 • 

The long-period Hamiltonian associated with each of these combinations of frequencies is of the 
form 


H*** (a, / 3 , t) - Ke j ~ 2 a ~ 3 cos + "3 ^3 + et ] - 

and the change of variables 


a 


* 

1 


a 


1 


0 ; 


6 t 

Z 3 ! + aT 
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leads to a Hamiltonian which is a constant of the motion 


ea* 

H ' = Kei /a x *a 2 a 3 + 6) a ^ 2 + w 3 /3 3 ) + — • 

The derivatives of a * , a 2 , and a 3 are linearly related and 


3 • * 


i.e 


•9 


a 3 = ^ ( a i* +C l) 


and 


i.e 


•9 


a 


2 



+ 


c 2 ) 


The Hamiltonian can now be written as 


H' 


Kei y y^a,* (a; + Cj ) ( a ; + C 2 ) cos (« t P* + * 2 P 2 + - 3 P 3 ) 


and, similar to the previous case, (t) is zero only when 


H 0 + — 

■V a i ( a i* +c i)( a r +c 2 ) = 


Ke J 


. V w 2 ~3 


In this case a* (t) is unbounded for any set of initial conditions such that a i (0) o-j (0) a 3 (0) ^ Oif 

^'K)L 


Ke^ yf Ccl a 
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where 


f ( a l* ) = ) / V« + C, ) « + C a ) . 

Thus, when e is sufficiently small, all three modes are excited and the pitch, roll, and yaw oscil- 
lations increase without bound. 

In a similar manner the combinations of frequencies shown as items 5, 7, 8, 18, and 27 of 
Table 5-1 are all of the form 

w 3 - 2co. ± j = 0 , i = 1,2 j = 0, 1, 2, 3 , 

and the long-period part of the Hamiltonian has the form 

H*** (a,/?, t) = Ke> a 1 cos (a> 3 P 3 - 2a>. /3 g + et ) , 

and, for j = 0, this is the internal resonance case of Reference 12. The analysis proceeds with 
new variables, 


p i - 


et 

co. 

i 


and a new Hamiltonian 


H' 


ea* 

Ke‘ a* r^r 3 cos (co $ P 3 - 2oj. / 3 .*) - • 


Then 


a* - - 2 co L Ke j a.* /aj cos (^ 3 /? 3 - 2co i ft*) 

and 


“3 " " 3 Kei a *^° r 3 c °s (« 3 ^ 3 -2a>. /?.*) . 
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Therefore 



0 


and 


a 


3 



= C . 


The new Hamiltonian can thus be written 


/ 

H' = Ke> a‘ yc~ ^ a* cos (o> 3 /3 3 “2a). J3') . 

The process that was demonstrated above applies also to this Hamiltonian. Both the phase plane 
and the geometric construction for finding values of a* for which a.* (t ) = 0 are shown in Ref- 
erence 12. In essence, the important features of this motion are that the sum of a.* and a 3 are 
constant as is a. (j / i) and, 


when 


_L 

2co. 


< Kei YC\ 


(5-5) 


there is an appreciable periodic interchange of energy between the pitch, a 3 , and roll-yaw a. (i = 1 
or 2) modes. 

In fact, it is possible to obtain an analytic solution for a.* (t). The first step is to eliminate 
the /3 r s from the equation for a* (t ) . The process is as follows: 


V(t) 


- 2o). Kei a.* j/c - a.* sin (a> 3 /3 3 ~ 2co. /3.*) 


or 


[i.*(t)] 2 - e 2 i a.* 2 ^C- 2 ^ a.*) [l- cos 2 (w 3 /3 3 - 2a). /J.)] . 


However, 


Kei °-i / c _ 2^ a i* c ° s ■ 2 "i ^i*) - 


ea. 

1 

2^7 
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and therefore 


or 


[“!*(*)] 


2 





1 

2 co. 


da.* 




/ 


2w. 


.* 3 + 


(^7 _CK2e2i ) a i* 2 + < +H 0 2 


(5-6) 


which is an elliptic integral of the first kind. This equation is applicable to any of the five cases 
listed above, whether they are internal or external resonances. 

The internal resonance for e = 0 and = 2^ was discovered numerically by Kane (Refer- 
ence 11). He considered a restricted form of three dimensional motion in which the roll and yaw 
angles were linearized and a circular orbit was assumed. Floquet theory was used to evaluate the 
stability of these equations at equal intervals over a parameter space (k x , k 2 ) where K x = - (l - r 2 )jr x 
and k 2 = (l - r r 2 . The initial pitch angle was the third parameter of his solutions. The pre- 
ceding analysis, when specialized to the case i = 1, j = 0, indicates that there will be a region 
about the line oj 3 = 2 o> x whose width is proportional to the initial pitch amplitude and in which there 
will be a large interchange of energy between the pitch and roll-yaw modes of oscillation. The re- 
sults of Kane's Floquet analysis for initial pitch angles of 5° and 10° are repeated here as Fig- 
ures 5-13 and 5-14 where these are taken directly from his Figures 8 and 9 respectively. The re- 
sults are replotted versus r x and r 2 in Figure 5-15 and as predicted in Reference 12 show very 
good agreement with the anticipated results. The scarcity of points at the two ends of the resonance 
line result from the non- uniformity of the test grid when it is transformed from (k x , K 2 ) to (r x , r 2 ) , 
and could be removed by evaluating the Floquet analysis for a finer grid. There are also three 
points of unstable motion on Figures 5-14 and 5-15 which are not explained by the above. These 
points may be the result of higher order resonance terms not included in this analysis. 

Kane also numerically integrated the exact equations of motion for several of his points of in- 
stability. His results show that these cases actually correspond to points at which the amplitude 
of the roll-yaw motion exhibits a slow but bounded oscillation as is apparent from the preceding 
analysis of the third order Hamiltonian. Breakwell and Pringle commented on this similarity and 
also demonstrated that the period of the variation in the amplitude of the roll-yaw oscillation ob- 
tained from Equation 5-6 was in good agreement with the result of Kane's numerical integration of 
the exact equations of motion. 

The combinations of frequencies shown as items 9, 17, 25, 26, and 31 in Table 5-1 are all of 
the form 

00 3 + 2co i " J = 0 ’ i - 1 , 2 j = 1 , • • • 4 ; 
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and the form of the Hamiltonian associated with these terms is 


H***(a,/3, t) - Ke J a. /a^ cos + 2w. “ j)t +cj 3 /? 3 + 2oj. /3.J . 

This general form resembles what resulted when the resonance co 3 ^ 2 a> x was evaluated in the Delp 
region in Reference 12. When the approach is generalized, the following results are apparent. The 
difference between the pitch variable a 3 and the appropriate roll-yaw variable a. is a constant, i.e., 


*3 <*) 


2^: a i (*) 


c . 


When the frequency mismatch e is large enough, bounded oscillations are possible for a 3 (t) and 
a. (t) . However, if | e /2 o> x Ke j | is sufficiently small, the motion of both a 3 (t) and a. (t) is com- 
pletely unbounded. Similarly the combinations of frequencies shown as items 6, 10, 11, 20, and 29 
are all of the form 

co + co. - co. ± k = 0 , i, j = 1 , i/j k~l,2, 3, 

1 3 

and they are all associated with terms in the Hamiltonian of form: 

H*** (a,/3, t) = Ke k ^ a 2 a 3 cos co 3 + oj. - co. ± k) t + a> 2 f$ 3 + a;. /3. - co. fi. J . 


In this case 


a . 

i 


co. 

i 



a . 

1 


C : and a 3 



C 


2 ■ 


The first of these two relations establishes that if the roll-yaw motion is initially small it will re- 
main small; the second shows that the pitch motion will also remain small. 

The frequency combinations shown as items 12 and 21 in Table 5-1 are both of the form 


co 3 - co x - co 2 + i ~ 0 , j - 1,2, 

and the associated Hamiltonian has the form 

H = Ke j cos [(“3 _Cl, l“" 2 + + 3^3 -a, l 
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The analysis for these two cases shows that 


1 


and 



C 


1 ’ 


a 3 a i 


Thus energy can be interchanged between the pitch a 3 and both of the roll yaw modes a l and a 2 . 

All three modes are bounded and there will be an appreciable interchange only if | e/co x Ke j | is suf- 
ficiently small. 

Application to DODGE Spacecraft 

The results of the preceding analysis find an interesting application in the DODGE satellite 
(Department of Defense Gravity Experiment) of the Johns Hopkins University Applied Physics Lab- 
oratory. This spacecraft was designed to test several combinations of mechanical configurations 
and damping methods. The flight experience of DODGE has displayed a number of unexplained 
anomalies in the actual versus predicted attitude motion (Reference 33). This analysis does not 
purport to explain these anomalies; however it does illustrate some interesting features. 

It is possible to vary the length of the inertia booms of the DODGE spacecraft and this feature 
has been exercised frequently during the first six months of the spacecrafts flight. The inertia 
parameters corresponding to two of the known configurations for DODGE are superimposed on 
Figures 5-16 and 5-17. At various times in the flight the configuration has been such as to map 
into other points in this same area. These two points straddle two of the lines of "unbounded" mo- 
tion which were developed in the preceding section. It is possible that for some portions 
of the DODGE Mission the inertia parameters have been such as to plot exactly on one of these 
lines. Actually, no spacecraft with similar extensible booms is exactly rigid; we cannot exactly 
equate the normal frequencies a> l9 co 2 , and w 3 with the telemetered boom lengths. Thus we cannot 
determine whether there have been times when the frequencies of the normal modes have exactly 
satisfied the criteria for resonance. 

Even so, there are several points that should be considered. Although, as a result of the un- 
certainty mentioned above, we cannot state conclusively that the otherwise unexplained large ampli- 
tude motions of the DODGE spacecraft are at least partly due to the resonance phenomena des- 
cribed herein, it is possible. These resonances depend on frequencies of excitation associated 
with terms of the third and fourth order in e. However, as mentioned in the previous section, the 
magnitude of the possible disturbance does not depend on the value of e. A numerical simulation of 
the spacecraft equations of motion that includes terms of third or fourth order in e from the devel- 
opments for the orbital parameters (as most do) cannot possibly model the effect of these resonances. 
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Figure 5-1 7— Resonance lines for large intermodal energy exchange 



The preferred method (this may sound like a paradox) of reducing any possible resonance ef- 
fects for DODGE is to shorten the main booms. This unfortunately tends to decrease the pitch and 
roll restoring torques, but it moves the parameter point (r lf r 2 ) toward a "safe" area of Fig- 
ures 5-16 and 5-17 (see below). The other alternative— lengthening the booms— leads to a more 
planar configuration; more possible resonances exist for this configuration than for any other, and 
it should be avoided if possible. 

Summary 

The motion of a nominally stable gravity -gradient- stabilized spacecraft in a nearly circular 
orbit is often described as consisting of small oscillations of its normal modes plus small forced 
oscillations due to the orbital eccentricity. This description is known to be inadequate if the in- 
plane (pitch) period is the same as the orbital period (oj 3 = 1); there is then a linear resonance be- 
tween the in-plane mode and the first term in the expansion for the non-uniform rate of rotation of 
the local vertical. In other respects, too, the situation is much more complex. Resonance rela- 
tionships exist whenever the normal frequencies are such that there are terms in the equations of 
motion for which 


ia j ± j <r^ 2 ± kc^‘ 3 - ±m, i , j , k , m , - 0 , 1 , ‘ * * , n . 

Thirty-two resonance combinations of frequencies for a single rigid body in a slightly eccentric 
orbit are shown in Table 5-1. These combinations were developed on the basis of a first-order 
solution to a set of equations expanded through the third order in the generalized coordinates and 
momenta and with terms through fourth order in e retained in the expansions for the orbital param- 
eters. Additional terms would undoubtedly appear if a higher-order solution were developed, or if 
more terms were included in the equations of motion or of the orbital expansions. 

Nineteen of the resonance combinations shown in Table 5-1 lead to unbounded increases in the 
amplitude of oscillation of one or more of the spacecrafts modes. Figure 5-16 plots the loci of 
inertias for which these resonances exist. Eight more of the combinations of frequencies lead to 
significant interchanges of energy between various modes of oscillation (e.g., a small pitch oscil- 
lation could change to a large roll-yaw oscillation and then back to a small pitch oscillation with a 
typical period for the interchange of twenty orbits). Figure 5-17 shows the loci of inertial combin- 
ations for which this occurs. The remaining five resonances have no significant effect on the motion. 

Several comments on the analysis of the preceding section are in order. As has been noted 
previously, the phrase "unbounded motion" means that the motion increases to a point for which 
higher-order terms become significant. A measure of the closeness of the frequency resonance 
of the form |e/Ke j w\ which is required to produce the various motions is developed for each situa- 
tion. The term K is a constant which is obtained by collecting terms in the expansion of the Ham- 
iltonian, as illustrated in the previous section for the case 2co 1 - 1 = 0. As such the value of k is 
dependent on the order to which the solution is formulated. Thus, these conditions for bounded 
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motion represent the relationship between e and the order of £ for which the resonance becomes 
important, and not the exact value of e, unless k includes higher-order terms. 

When the inertia parameters of a spacecraft lie exactly on any of the lines of Figures 5-16 or 
5-17, the motion will be resonant no matter how small e is. The rate of change of the appropriate 
mode will depend on e, but the final value will not. This was also demonstrated in Chapter 2, where 
each of the unstable regions exists for an arbitrarily small e. The dependency of the rate of di- 
vergence on e is indicated for the planar case by Figures 2-6 through 2-8. 
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CHAPTER 6 


CONCLUSION 


Summary 

The analyses in the preceding chapters have demonstrated a number of differences between the 
motion of a rigid gravity- gradient- stabilized spacecraft as determined by a linearized analysis for 
a circular orbit and as determined by a more complete analysis for the same spacecraft in a slightly 
eccentric orbit. Prior investigation, both analytical and numerical, had previously indicated the 
importance of considering the various effects studied herein. The major contribution of this dis- 
sertation is the specific demonstration— in a manner suitable for use in spacecraft design— of the 
effect and location of both parametric and nonlinear resonances. 

There are a large number of resonances that can affect the motion of a gravity- gradient- 
stabilized spacecraft. Several possible resonances occur for inertia parameters in the portion of 
parameter space at present considered most desirable for spacecraft design. Many of the res- 
onances only occur for eccentric orbits, and several are the result of terms which are third- or 
fourth-order in eccentricity. The rate at which a particular resonance affects spacecraft motion 
is a function of the magnitude of the eccentricity as is the width of the band in parameter space for 
which the resonance occurs. However, in any given resonant situation the peak magnitude of the 
oscillation does not depend on the value of the eccentricity, although for the energy-interchange 
types of resonance it does depend on the initial conditions. 

Recommendations 

The existence of resonances, such as those studied in the preceding chapters, has consequences 
that should be considered in the design of any gravity-gradient-stabilized spacecraft. In the con- 
ceptual design stage, configurations which deliberately utilize an interchange type of resonance to 
transfer energy from one mode to another to enhance damping (e.g. Pringle (Reference 34)) should 
be considered. Once any configuration has been selected, the parameter selection process should 
specifically include consideration of all the types of resonance. 

In general, much of the analysis of a given spacecraft is based on the results of a digital simu- 
lation of the equations of motion. These simulations frequently retain the nonlinear character of the 
equations of motion as is desirable. However, if the expansions for the radius vector and true anom- 
aly are limited to first-order terms in eccentricity, the numerical simulations cannot demonstrate 
the effect of many of the pertinent resonances. Terms through the fourth order in eccentricity should 
therefore be included in the simulations even though in nonresonant situations they have little effect. 
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Extensions 


The analysis of this dissertation could be extended in several directions: it could include 
higher-order terms in the Hamiltonian, or continue the method of averaging to second order. How- 
ever, the return from either of these steps is not likely to justify the effort. 

The analysis could also be extended to cover several types of spacecraft presently being.con- 
sidered. The equations could be quickly modified to include the effect of a constant- speed inertia 
wheel as has been suggested for improved yaw response in a number of spacecraft. Or the effects 
of damping on the motion of an essentially rigid spacecraft (achieved through use of hysteresis rods 
or a greaseball type of damper) could be considered. The order of the system could be extended 
to include two- or three-body spacecraft, thus including configurations such as the Vertistat. 

The analysis of the fourth and fifth chapters follows a formal procedure which is applicable 
to a large class of mechanical systems whose motion consists of "small oscillations.” The Ham- 
iltonian for any such system can be written as the sum of terms of increasing order in the coordi- 
nates and momenta (i.e. H r h 0 + H x + * * • + H n ) . It may be possible to determine generalized criteria 
for parametric and nonlinear resonances from the form of the third and higher order terms in the 
Hamiltonian without requiring a complete analysis such as the preceding. For example if in the 
second order Hamiltonian there exists a coordinate pair q £ , p. such that there are no terms in 
q i Qj 9 Qi Pj or Pi Pj f° r * ^ i 9 then the linear equations for the q. and p. are independent of the 
remaining variables. There may be similar criteria for the third order Hamiltonian which de- 
termine the existence or absence of the various resonances. 

Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, November 5, 1969 
039-01-01-01-51 
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Appendix A 


Basic Equations of Motion 


The motion of a rigid spacecraft in an inverse square force field can be completely specified 
by six coordinates. The position of the center of mass can be specified by three coordinates r, p, 
v and the attitude of the body about its center of mass by the Euler angles <f>, 0, and \p. The equa- 
tions of motion for this six -degree -of -freedom problem can be obtained directly from the potential 
and kinetic energies via the Lagrangian formalism. 

The kinetic energy of the spacecraft is 

T T m s [r 2 +( r P) 2 + (rcospv) 2 ] +y [li^i 2 


+ 1 2^ + h^ 2 ] 


where 

0 1 = 0-(^>+v)sin(P+p), 

Q 2 r (b + p) cos Kp + (cp + v) cos ( 0 + /9) s in 0 , 

and 


fi 3 - - (@ + p) simp + (<0 + v) cos (P + p) cos \p 


The potential energy of the spacecraft is 


£l + ^ 3 ^ 2 3 

r 2r 2 + 2r 2 


|^(l 2 - Ij) (cos cp cos 0) 2 


+ 




(cos <j> sin 0 cost/f + sin <£ s in 
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The three Lagrangian equations for the coordinates describing the motion of the c enter -of -mass 
of the spacecraft are 


m s ir* - m s r p 2 - m s r cos 2 pv 2 — ~ i..* s 


1 1 + *3 * 2 


m + 3 

r ■* | s 2r 2 


^ 2 ^(1 2 ” * i) ( cos 0cos 0) 2 + (l 2 - I 3 ) (cos sin 0 cos i/>+sin0sim//) 2 S - 0 


f 


dt -< m_ r 2 cos 2 


p v - 1 1 - (<£ + v ) sin (6? + p)J sin (0 + p) 

+ p) cos ip + (<p + v ) cos + p) sin cos (0 + p)sini/> 

+ I 3 (6' + p)sin0 + (0 + v) cos (0 + p) cos i//J cos (£ + p) cos ^ - 0 


+ In 


and 


(m s r 2 p + I 2 co 2 cos 0 - 1 3 ^3 s i n 0 ) + m s r 2 v 2 cos p sinp 

+ I A a) 1 (<£ + v) cos (0 + p) + I 2 C0 2 + v ) s in \p sin (0 + p) 


+ 1 3 (<£ + v ) cos P cos + p) - 0 . 

The three preceding equations all contain terms of sharply contrasting magnitudes since the in- 
ertias I 1? l 2 , and I 3 are all of the order m s ^ 2 where l is the linear dimension of the spacecraft. 

The ratio (Z/r ) 2 is on the order of 10" 12 for a gravity gradient spacecraft, thus the equations for 
the orbital coordinates can be simplified to the form 


V - r p 2 


r 


[< 


cos p V ) 


2 




0 



0 


dt [ f2 ^] - 0 

and these are the well known equations for the motion of a point mass subject to an inverse square 
gravitational force. Actually Beletskii demonstrated in Reference 9 that the relativistic correction 
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due to the velocity of a typical gravity stabilized spacecraft is several orders larger then the 
terms deleted above, although it too can be safely ignored. The coordinate system can be chosen 
such that p = p - o in the above equation leaving the two familiar equations 

p, 

r - r v 2 + — ~ ~ 0, 
r 2. 


and 



= 0. 


The equations for the remaining coordinates are 

1 1 ^0 - (0 + v) sin 0 - (0 + v) 6 cos + (l 3 “ I 2 ) |(0 + v) cos 0 cos 20 


+ \ (0 + v) 2 cos 2 @ + sin 2 24 > + ~2 @ 2 sin 2 0 


3/x 


~3~ {^3 ~ ^2) (sIh 2 0 sin 20 + cos 2 0 sin 2 6 sin 20 


sin 20 sin 0 cos 20) - 0 


(l 2 cos 2 0 + 1 3 sin 2 0)0 + (l 3 - I 2 ) 


9 \p sin 20 — 2~ (0 + v) cos 0 sin 2 0 


+ ~2 (0 + v) 0 s in 6 s in 20 ~ (0 + 


v) 0 cos 6 cos 2i/J - (0 + v) 2 sin 2 d 

- (0 + v) 0cos + ~2 (l 2 " I 3 ) (0 + v) 0 sin 20 sin 20 


1 

+ 2 


(l 2 sin 2 0 + I 3 cos 2 ^) (0 + v) 2 sin 2 0 


+ 



1 t ) cos 2 0 sin 20 - (i 2 - 1 3 ) (cos 2 0 s in 2(9 cos 2 0 


1 

+ 2 


sin 2^ cos 0 sin 20 


>] 


0 
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and 


(i 3 cos 2 0 + I 2 sin 2 0) 


£(0 + v) cos 2 @ - (4> +v) 9 sin 2#j 

1 £(0 + v) sin 2 9 + (0 + v) 6 sin 29 - 0 sin 9 -00 cos 


+ I 


+ ~2 (l 2 — I 3 ) cos @ sin 20 “ 9 2 sin 9 sin 20 
+ 2^0 cos 9 cos 20 + (0 + v ) 0 cos 2 <9 sin 20^ 


3fJ. 


^I 2 ~ I ^ sin 20 cos 2 9 - ( 1 2 — 1 3 ) (cos 20 sin 6 sin 20 


+ sin 20 sin 2 0 - sin 20 sin 2 9 cos 2 0)j - 0 . 
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Appendix B 


Auxiliary Formulas Relating to Orbital Mechanics 


A collection of formulas which are applicable to the motion of a point mass in an inverse 
square gravitational field are given below. 

Force of Gravity 


-F 


yMm 

“ 72 ~ 



fj . m 


r 2 


Motion in an Elliptical Orbit 


p 

T + e cos v 


Useful Derivatives 


p = b [l - e 2 ] 1/2 = 


r 


2 


dv 

dt 


K 


KT 



2t t ab 


K _ 277 

ab T 


K 2 a _ i± 
a 3 b 2 " I 3 


dv 

dt 


K 

2 


K( 1 + e cos v) 2 


d dv 
dv dt 


dv _ ~2eK sin v 
dv ~ rp 
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Change of Independent Variable 


d[ ] 

d [ ] 

dv 

dt 

dv 

dt 

d 2 [ ] 

d 2 [ ] 

/dv 

dt 2 

dv 2 

\dt, 


K 2 d 2 [ 

J _ 


r 4 dv 2 


JK d[ ] 
r 2 dv 


2 d [ ] dv dv 
+ dv dv dt 


2eK 2 sinv d [ ] 
pr 3 dt 


At( 1 + e cos v ) d 2 [ ] 
r 3 dv 2 


M 2e s in v d [ ] 
^3 dv 
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Appendix C 


Development of Asymptotic Expansion Formulas to the Third Order 

The method of application of asymptotic expansions to problems of nonlinear oscillations 
which is used herein was developed by the Russian authors N. M. Krylov and N. M. Bogoliubov 
in Reference 26. A more thorough explanation of the method is given by N. M. Bogoliubov and 
Y. A. Mitropolsky in Reference 27, and readers interested in the theory underlying the develop- 
ment of the formulas which follow should consult the latter reference. A considerable simpli- 
fication in the development results from the restricted form of the disturbing function for this 
problem. 

The intent of this approach is to find a solution to the differential equation 

d 2 X 

— + ^ X = £ f ( t , x) (C-l) 


in the form 


m~ 1 

x - a cos ■/' + ^ * e n u p (a, >fi, t ) 

n = 1 


(C-2) 


where a and 0 are variables obtained from the equations 


m 

da V -1 

dF = 2- e"A n 


(a) 


(C-3) 


and 


d0 

dt 


a, + 


m 

L. '" B " 


(a) 


(C-4) 


The perturbing function which is needed for this application is of the form 


ef( 


4 

t , X) = £ e n [ f nl (0 + f n2 (t)x] . 


(C-5) 
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Formally the problem becomes one of finding functions A n ( a ), B n (a), u n (0, 0, t) n = 1, • • ■ , m 
such that Equation 02, with a and 0 defined by Equations 03 and 04, satisfies Equation Ol to 
an order e m+1 . This is accomplished by differentiating Equation 02 twice with the aid of Equa- 
tion 03 and 04, and substituting the result into the left hand side of Equation Ol. Then the 
right hand side of Equation Ol is also expanded by substituting Equation 02 for x in f(t, x). 
Terms of equal magnitude in e are then equated to give a series of differential equations from 
which the desired functions are obtained. The resulting equations are all of the form 


d\p 


2 + 2a d0dt + 


= F m + 2cc A m si.rn// + 2aa ; B m cos 0 (C-6) 


where 


F 1 ” f u( t ) + f 12( t ) a COS 'P 


F 2 = f 2 i (t) + f 22 (t) acos</> + f 12 (t)uj 


( M 

+ yaBj 2 - Aj Aa J cos + aAJ da 


+ 2Aj I sin 0 


d 2 u x d 2 u x d 2 Uj ^ 2 

2 ^>B 1 0^2 ~ dadt ~ ^ B i d\pdt ~ \ da 


F 3 = f 31 (t) + f 32 aCOS ^ + f 22 C 4 ) 11 ! + f 12 u 2 


/ dA 2 dA t \ 

y2a Bj B 2 ~ Aj da ~ ^ 2 da J 


dB 2 dBj 

cos 4 , + laA, -gr + aA 2 ~ A ^ 


+ 2Aj B 2 + 2A 2 B, ) sin ip - A 


\ / dAj ^u 1 dB 1 du^ 

ij s i n 0 A i ^ da (9a + da <90^ 


<9 2 u d 2 u x <9 2 Uj a 2 u 1 

- A x 2 “ (B x 2 + 2 ojB 2 ) ^2 ~ ( 2A i B i +2a ' A 2 ) d a d0 ” 2A 2 da <9 1 


d 2 d 2 u 2 

2B 2 d0dt dadt ~ 2a; A i dad0 ” 2cd B i 3^2 ” 2B 1 d0dt 


d 2 u 2 d 2 u 2 d 2 u 2 
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and 


F 4 - f 41 <*) + f 42 (t) 3COS0 + f 32 (t) Uj + f 22 (t)u 2 + f 12 (t) Uj 


aB 2 2 -A — ‘ — ' 


( 2aB i B 3 + 

^ aA l da + 

/ dA^ dA t \ <9u 1 / dB 2 

\ A l da + A 2 da/ <9 a ” l A i da 


dAA 

‘i ”dIT " A 2 "diT A 3 da / cos,/ ' 


dB 2 dB, 

aA 2 "da” + aA 3 di” + 2Aj B 3 + 


2A 2 B 2 +A 3 B i) 


sin 0 


HR \ 


d 2 u 

2A > a 2 7 ^ 


d 2 - 


- 2 B 


3 dt dip 


d2 u i 


(2B» B 2 

u ^ 

+ 2 “ B 3) 3^ 

- 2 (A x B 2 + A 2 B 

d 2 u x 

dA t 5 u 2 

dB l 3U 2 

; dt 2 

~ A i da d a 

A i da (90 



A ! 2 ^” ( B ! 2 + 2 -B 2 ) 


d 2 u„ 




( 2 Ai B , + 2 “^)-^ 


2A 


2 <9a (9t 


- 2B 


2 (90 < 5 1 


d 2 u 3 <9 2 U_ 

- 2co B, . — • - 2 w, A, 


(90 


"l n i (9 a (90 


< 9 2 u 3 (9 2 u 3 

2A i da<9t " 2B i d0dt * 

The first equation is solved for u 19 A x , and B t as follows. Assume that u x and F t can be written as 
double Fourier series, 


u ! ( a > 0» t) 


Z2> ( a) exp [i(nt + m0)J , 


and 


Fi (a, 0, t) = 

B nm ( 3 ) ex P [i( nt + m 0 )] > 


n m 


117 


equate the coefficients of equal harmonics on both sides of the first of Equations C-6 to obtain 


F 


i 

n m 


( n + moj) 2 


a? - (n+mo;) 2 / 0 


and require that 


A x sin i p + 2awBj cos \)j + 


exp [i(nt + 


0 


for all n, m such that to 2 - (n +m cS) 2 = 0. This development is only valid for the M nonresonance case” 
which is characterized by a requirement that 


CO 



for integer values of p and q (other than p - q - l) such that 


_p _ n 
q 1 - m 

for values of n and m that occur in the expansion of F 1 as a Fourier series. (A similar set of 
formulas, which are valid for the "resonance case," are developed later in this appendix.) 
The only terms in F^ for which oj 2 - (n + iW) 2 m 0, subject to the above restriction, occur for 
n = 0, m = ±1. This part of F x can be written as 


F 1 (a,i//, t) - hj(a)sirn/i + h 2 (a)cosi// 


and Aj (a) and B x (a) are 


hj (a) 

M a ) - 2u 


and 


h 2 (a) 

B i ' 
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The values of A x , b x , and u l obtained from the above process can be used to expand F 2 (a, 0, t ) as 
a Fourier series from which A 2 , B 2 , and u 2 can be found and this process can be continued to any 
desired degree of approximation. 

The first order solution to Equation C-l is expressed as 

x = a co s \p 


where 


da 

dt ~ eA i o) 


and 


ddj 

d F “ + eB i (a) - 


the m th order solution is 


X 


m- 1 



n = 1 


e n u n (a, 0, t) 


where 


da V -1 

dF = 2^ e n A n (a) 

n = 1 

and 


d \p 

dt 


OJ + 2^ B n ( 3 ) • 

n = 1 


Several changes are made in the above process for the '’resonance case ,T which occur when 
p/q. In this case A is defined by 


eA = 
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Equation C-l is rewritten in the form 


"dr + (■?■) x = e [ f ( 4 - x)-ax] 


dt 2 V q 

and there is a change of variables from 4> to 6 where 


* -ft + 0 


The form of the solution becomes 


x - a cos 


t + <?) + «" U n O' t) 


where a and 6 are obtained from 


da 

dt 


m 


(a, 0) 


and 


d 6 _ 

dt * L 


e n B (a, 6 ) . 


The disturbing function is still in the form of Equation C-5 and the set of Equations C-6 is re- 
placed by 


d 2 u 
dt 2 


1 + m U 1 = F 1 ( a > ^ O + 2 |" A 1 sin ^ 


( 2 ? b i - a ) 


a cos \p , 


d 2 U 

d 


— + /-P V 

t 2 W 2 


- F, (a, </>, t) + 


r 

|2Aj Bj + a ^ a Aj 


1^1 P , 

+ a B i + 2 ~q A 2 | Slr "^ + 




dA 1 
da A i 


dA 


W B j + a Bj 2 + 2a B 2 cos i p 
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and 


?? u jl . (P \ 2 
at 2 U) U3 


F 3 (a, t//, t) + 2A 1 B 2 +2A 2 B 1 + a'^-A 2 


f [ 2A ! 

3B 2 p I 

IF B ! + 2 q a 3 J 31 

r <3Aj 5A 2 3 Aj 

|2aBj B 2 5a A 2 _ 5a A i 55 


5B 2 5Bj 

+ a “5T Aj + a -Q 0 - B 2 + a 


sin t// 


^2 p 

" 55 B l + 2 q aB 3 


where 

F x ( a , •>, t ) = f 11 (t) + f 12 (t)acos0 > 


F 2 (a, i//, t) - f 21 (t) + f 22 (t) a cos^; 


5 2 Uj 5 2 u j 

+ [f 12 (t) - A] u t (a, (/., t) - 2A, JJJt - 2B, jg^t 


F 3 (a, i p, t) = f 3 j ( t ) + f 32 (t) a cos i p + f 22 C t )u 1 (a, i p, t) 

fd k 1 5A, \ 5u 2 

+ [f 12 (t)-AJ u 2 (a, t/>, t) - \^ a Aj + d0 B l J dg 

fdB 1 9B t \ 5Uj 5 2 u t 

“ y5a~ a j + 'Je B, J -Jg- - A x 2 

5 2 Uj 5 2 u t 5 2 u 2 

2Aj Bj g a g@ 2Aj 2A t 


5 2 Uj 5 2 Uj 

~ B i 2 36 2 ~ 2B 2 55 5 1 ” 


COS 0 


<9 2 u 2 
2B i JdHt 
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\ 


The rest of the procedure is quite similar to that for the nonresonant case. 


u 1 is 

nm 


However this time 


F 1 


whenever 

(O’ * 

and u n J m = 0 otherwise. The previous inequality could also be written as 

nq+(m±l)p / 0 

and there are a number of values of n, m for which it is not satisfied (e.g., when p = 1, q = 2 it fails 
for n = 1, m = -1). The terms in F x (a, 0, t) for which the inequality fails can be written in the form 

hj (a, 0, t) sin 0 + h 2 (a, 0, t) cos 0 



and by equating the harmonics in sin 0 and cos 0 , 


\ (a, 0, t) 


Tp“ h i ( a ’ & 


t) 


and 


B 1 (a, 0, t) 



h 2 (a, 0, t) 
a 


As in the nonresonant case this process can be continued to obtain A n , B n , and u n to any desired order. 
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